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PREFACE 
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CHAPTER  1.   EXECUTIVE  SUMMARY 


The  objective  of  this  part  of  the  FAST/revenue  service  correlation  study 
was  the  development  and  application  of  the  methodology  for  prediction  of  the 
geotechnical  aspects  of  track  performance.   Track  performance  is  defined  in 
terms  of  vertical  track  settlement.   The  following  field  sites  were  evalua- 
ted:  1)  in-service  revenue  track  having  concrete  crossties  at  Streator 
(Leeds),  Illinois,  Richmond  (Lorraine),  Virginia,  and  Aberdeen,  Maryland,  2)  a 
revenue  wood  tie  track  section  at  Streator  (Leeds),  Illinois,  and  3)  wood  and 
concrete  tie  sections  at  the  Department  of  Transportation  Facility  for  Accele- 
rated Service  Testing  (FAST)  test  track  in  Pueblo,  Colorado.   The  Streator  site 
is   owned  by  the  Atchison,  Topeka  and  Santa  Fe  Railroad.   The  Richmond  section  is 
owned  by  the  Chessie  System.   The  Aberdeen  test  section  is  owned  by  Amtrak. 

A  variety  of  field  tests  were  conducted  at  these  sites  to  assess  the 
physical  states  of  the  ballast  materials  as  they  existed  prior  to  a  track 
maintenance  (raising  and  tamping)  operation.   The  ballast  physical  states 
were  identified  using  the  plate  load  test,  ballast  density  test, and  lateral 
tie  push  test.   These  physical  state  tests  were  repeated  after  the  surfacing 
operation  to  determine  the  changes  in  the  ballast  physical  properties  caused 
by  the  maintenance.   One  of  the  primary  effects  of  track  maintenance  was  a 
loosening  and  disturbance  of  the  traffic-compacted  and  stabilized  ballast 
layer,  both  in  the  cribs  and  under  the  ties.   The  physical  states  of  the  bal- 
last found  immediately  after  the  maintenance  were  approximately  independent  of 
ballast  type  and  initial  condition. 

Ballast  materials  were  recovered  from  the  sites  and  subjected  to  repeated 
load  triaxial  testing  to  investigate  the  elastic  and  inelastic  behavior.   A 
granite  ballast  from  FAST  was  also  tested  in  a  variety  of  repeated  load  condi- 
tions, including  a  mix  of  load  levels. 

The  elastic  behavior  was  defined  in  terms  of  the  resilient  modulus.   The 
resilient  moduli  values  from  extension  tests  were  smaller  than  the  moduli 
values  from  compression  tests  under  similar  stress  conditions.   When  the  ver- 
tical stresses  were  cycled  between  extension  and  compression  states,  the  total 
resilient  moduli  fell  between  the  individual  extension  and  compressive  values. 
The  resilient  behavior  of  the  field  site  ballasts  was  not  significantly  depen- 
dent upon  ballast  type. 

The  inelastic  behavior  of  the  ballast  was  found  to  be  dependent  upon  the 
first  cycle  strain,  Z. ,  and  a  material  constant,  C.   This  constant  was  inde- 
pendent of  stress      conditions  for  a  particular  ballast,  but  was  not  the 
same  for  all  ballasts.   The  total  permanent  or  inelastic  strain  developed  in 
a  ballast  sample  with  a  mix  of  loads  was  controlled  by  the  maximum  past  load 
and  was  independent  of  the  loading  sequence.   Subsequent  applications  of  re- 
duced load  as  well  as  partial  unloadings  resulted  in  negligible  increases  in 
permanent  strain.   This  means  that  the  largest  dynamic  wheel  loads  cause  most 
of  the  track  settlement.   A  strain  superposition  method  was  developed  to  ac- 
count for  the  deformations  that  result  from  repeated  applications  of  the 
mixed  loading  conditions  present  at  the  field  sites. 

The  inelastic  behavior  of  ballast  was  also  investigated  using  a  box-type 
testing  device,  intended  to  simulate  field  behavior.   The  trends  in  ballast 
deformation  from  these  box  tests  were  similar  to  those  observed  from  the  tri- 
axial tests.   The  box  tests  also  showed  that  relatively  large  horizontal 
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residual  stresses  can  develop  in  the  ballast,  both  in  the  crib  and  under  the 
ties,  as  a  result  of  cyclic  loadings. 

Measurements  of  vertical  track  modulus  at  the  revenue  sites  did  not  indi- 
cate that  significant  changes  in  the  magnitude  of  track  modulus  occurred  due 
to  the  maintenance.   The  major  factors  contributing  to  the  magnitude  of  track 
modulus  were  ballast  depth  and  subgrade  stiffness.   Also,  the  field  measure- 
ments of  track  modulus  showed  that  the  uniformity  of  track  support  was  not 
improved  as  a  result  of  the  surfacings. 

Among  the  important  differences  between  the  sites  was  the  train  loading 
imposed  on  the  track  systems.   The  dynamic  wheel  loads  at  each  site  were 
measured  by  BCL  to  characterize  the  mixed  loading  conditions.   The  static 
wheel  load  distributions  for  each  site  were  obtained  from  consist  data.   Al- 
though the  average  static  and  dynamic  wheel  loads  were  about  equal  for  a  par- 
ticular site,  the  static  wheel  load  distributions  for  the  revenue  sites  were 
quite  different  from  the  dynamic  distributions.   One  out  of  one  thousand  (0.1 
percent  exceedance  level)  dynamic  wheel  loads  at  Aberdeen  was  measured  at  75 
kips.   At  the  0.1  percent  exceedance  level  at  Leeds,  Lorraine, and  FAST,  the 
dynamic  wheel  loads  were  about  46  kips.   These  very  high  dynamic  loads  at 
Aberdeen  were  a  result  of  a  combined  effect  of  wheel  tread  irregularities  and 
the  high  speeds  of  passenger  trains.   The  loading  environment  at  FAST  was  much 
more  uniform  than  at  the  revenue  sites,  because  of  the  controlled  nature  of  the  FAST 
experiments.   The  average  static  and  dynamic  wheel  loads  at  FAST,  however, 
were  about  1.5  to  2  times  larger  than  the  average  values  measured  at  the  reve- 
nue sites. 

Survey  measurements  were  made  at  each  of  the  field  sites  by  BCL  to  deter- 
mine cumulative  vertical  track  settlements.   The  survey  data  from  Aberdeen  and 
FAST  were  supplemented  by  inductance  coil  strain  gauges  embedded  in  the  bal- 
last layer.   The  strain  gauge  instrumentation  of  the  ballast  layers  lead  to 
the  conclusion  that  the  zone  of  ballast  disturbed  by  the  maintenance  operation 
was  the  primary  source  of  track  settlement,  and  that  over  the  relatively  short 
experimental  time  period,  the  subgrade  contributions  to  track  deformation 
could  be  neglected  for  these  sites. 

The  settlements  at  the  Lorraine  and  Aberdeen  sites  continued  to  develop 
at  what  appeared  to  be  a  steady  rate  throughout  the  period  of  the  measure- 
ments, after  the  initial  traffic  following  maintenance.   Although  the  total 
vertical  settlements  at  all  of  the  revenue  sites  continued  to  increase  as  tonnage 
accumulated,  the  roughness  of  the  track  profiles  did  not  increase. 

The  ballast  deformations  at  FAST  were  found  to  be  similar  for  both  wood 
and  concrete  tie  sections.   This  was  true  for  the  ballast  measurements  obtain- 
ed during  the  initial  FAST  experiments  and  for  the  ballast  deformations  fol- 
lowing the  1979  rebuild. 

A  three-dimensional,  nonlinear,  elastic,  multilayer  model,  GEOTRACK,  was 
used  to  predict  track  deflections  and  track  moduli  for  the  revenue  sites.   The 
model  predictions  of  vertical  track  deflection  agreed  well  with  field  meas- 
urements, although  the  model  tended  to  underestimate  the  field  values. 

The  methodology  to  predict  the  inelastic  or  permanent  settlement  of  the 
field  track  sections  used  GEOTRACK  nonlinear  stress  analyses  of  the  field  site 
track  sections  with  stress-dependent  ballast  properties,  representative  sub- 
grade  moduli,  and  multiple-axle  combinations  of  the  measured  dynamic  wheel 
loads.   The  predicted  stresses  from  the  computer  analyses  for  each  particular 
field  site  and  loading  condition  were  used,  along  with  the  permanent  deforma- 
tion behavior  of  the  ballasts  from  laboratory  tests  to  estimate  the  resulting 
ballast  strains. 
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The  general  magnitudes  of  the  predicted  deformations  agreed  with  the 
measured  values.   The  shapes  of  the  deformation  versus  million  gross  tons 
predicted  results  did  not  agree  as  well.   The  predicted  rate  of  deformation 
accumulation  for  the  revenue  sites  was  less  than  the  survey  measurements 
indicated,  particularly  for  the  Aberdeen  and  Lorraine  sites.   The  settlement 
predictions  made  for  FAST  Section  22  agreed  best  with  the  field  measurements. 
This  was  probably  a  result  of  the  more  controlled  track  loading  and  test 
conditions  present  at  FAST. 

Uncertainties  in  both  the  first  cycle  strains  and  the  representative 
stress  paths  could  account  for  some  of  the  differences  between  measured  and 
predicted  settlements.   Deviations  between  measured  and  predicted  track  set- 
tlements could  have  resulted  from  other  factors  not  considered  in  the  method- 
ology.  One  of  the  most  important  factors  not  included  is  thought  to  be  set- 
tlement induced  by  vibrations  transmitted  through  the  track  superstructure  to 
the  ballast.   Other  factors  that  can  have  varying  degrees  of  influence  on 
track  settlement  and  performance  include  ballast  degradation,  environmental 
effects  and  subgrade  responses.   Although  these  last  three  factors  did  not 
appear  significant  for  the  monitoring  periods  at  the  field  sites  investiga- 
ted during  this  study,  they  could  be  significant  at  other  sites. 
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CHAPTER  2.   INTRODUCTION 


In  order  to  improve  current  track  design  methods  and  quantify  track  main- 
tenance life,  research  is  necessary  to  understand  the  behavior  of  the  track 
foundation  under  dynamic  traffic  loading.   The  repetitive  loading  and  unload- 
ing of  track  structures  due  to  the  passage  of  train  traffic  causes  settlement 
in  the  supporting  ballast  and  subgrade.   As  the  traffic  history  accumulates, 
these  deformations  continue  to  increase  until  the  track  has  deteriorated  to  a 
state  at  which  it  can  no  longer  perform  at  a  desired  service  level.   At  this 
point,  the  track  is  surfaced  (raised  and  tamped)  to  reestablish  the  necessary 
profile,  and  the  deformation  process  starts  again. 

The  forcus  of  this  report  is  on  the  prediction  of  track  performance  at 
four  field  sites.   Track  performance  is  assessed  in  terms  of  vertical  track 
settlement  and  the  factors  that  contribute  to  settlement.   The  sites  consist 
of  concrete  tie  revenue  service  track  at  Leeds,  near  Streator,  Illinois;  Lor- 
raine, near  Richmond,  Virginia;  a  site  south  of  Aberdeen,  Maryland;  and  wood 
tie  track  at  Leeds.   Wood  and  concrete  tie  track  at  FAST  is  included  for  com- 
parison. 

Differences  in  track  structural  characteristics  and  track  foundation 
among  the  sites  were  identified,   A  variety  of  field  tests  were  conducted  to 
assess  the  physical  states  of  the  ballast  materials  as  they  existed  prior  to 
track  maintenance.   These  physical  state  tests  were  repeated  after  a  surfac- 
ing operation  to  determine  the  changes  in  the  ballast  properties  that  occur- 
red due  to  track  maintenance.   To  complete  the  site  characterization,  stan- 
dard penetration  and  Dutch  cone  penetration  tests  were  used  to  investigate 
the  subgrade  support  characteristics  of  the  revenue  sites.   These  tests  were 
done  using  a  drill  rig  specially  equipped  with  rail  wheels  so  that  the  borings 
could  be  taken  directly  under  the  track  structure.   The  results  of  these  sub- 
surface explorations  were  used  to  help  identify  variations  in  strength  and 
composition  of  the  foundation  layers  over  the  depths  influenced  by  train  load- 
ings.  Moisture  changes  in  the  top  of  subgrade  were  also  monitored  and  corre- 
lated with  weather  conditions.   A  summary  evaluation  of  these  test  results 
has  been  given  in  a  companion  Volume  I  [1]*  to  this  report.   The  results  of 
these  tests  will  be  explained  further  where  necessary  to  support  conclusions 
made  within  this  report. 

Among  the  important  differences  between  the  revenue  sites  are  the  train 
loadings  imposed  on  the  track  systems.   To  determine  these,  the  dynamic  wheel 
load  distributions  were  measured  by  the  prime  contractor,  Battelle-Columbus 
Laboratories  (BCL) ,  for  representative  traffic  at  each  of  the  sites.   Since 
each  site  was  subjected  to  a  significant  range  of  wheel  loads  rather  than 
constant  repetition  of  the  same  magnitude  loading,  a  method  was  developed 
that  can  be  used  to  account  for  the  real  mix  of  wheel  loads. 

The  ballast  and  subgrade  soils  from  each  of  the  field  sites  were  tested 
to  determine  both  the  resilient  and  inelastic  behavior  of  these  materials 
under  representative  loading  conditions.   Both  cyclic  and  repeated  load  tri- 
axial  tests  were  used  to  determine  the  resilient  modulus  of  the  ballast  ma- 
terial under  full  unloading,  partial  unloading,  and  shear  stress  reversal 
conditions.   Undisturbed  subgrade  samples  were  also  subjected  to  repeated 
load  tests  under  anticipated  field  stress  conditions  to  determine  the 


^Numbers  in  brackets  [  ]  denote  references. 

12 


resilient  modulus.   The  final  estimates  of  the  elastic  properties  of  the 
track  substructure  were  based  on  the  results  of  the  physical  state,  triaxial, 
standard  penetration  and  cone  penetration  tests. 

The  development  of  an  analytical  model  that  incorporates  the  fundamental 
track  structural  components  and  foundation  characteristics  is  a  necessary  tool 
for  assessing  track  performance.   The  model  used  in  this  study,  GEOTRACK,  is 
a  three-dimensional,  nonlinear,  elastic,  multi-layer  track  analysis  program 
that  emphasizes  the  geotechnical  aspects  of  track.   The  GEOTRACK  model  is  con- 
ceptually similar  to  the  MULTA  model  developed  by  BCL,  as  explained  in  Ref . 
[40]. 

Measurements  of  vertical  track  modulus  were  taken  before  and  after  main- 
tenance at  the  field  sites  by  BCL  during  the  initial  site  visits.   The  track 
structural  and  foundation  properties  determined  for  the  field  sites  were  used 
with  the  GEOTRACK  model  to  predict  values  of  track  modulus.   The  comparisons 
between  measured  and  predicted  elastic  track  moduli  served  as  an  additional 
validation  of  the  model's  elastic  capabilities. 

The  prediction  of  the  vertical  track  settlement  at  each  of  the  field 
sites  was  based  upon  inelastic  behavior  of  the  foundation  soils,  as  measured 
in  the  laboratory  under  representative  loading  conditions.   The  inelastic  be- 
havior of  the  foundation  soils  has  been  investigated  using  repeated  load  tri- 
axial tests  and  a  ballast  box  device  developed  for  this  program  for  simulating 
field  conditions.   A  method  was  developed  that  can  be  used  to  account  for  the 
deformations  that  occur  from  repeated  applications  of  the  mixed  loading  condi- 
tions present  at  the  field  sites.   Predicted  values  of  track  settlement  were 
then  compared  with  the  track  settlements  measured  at  the  field  sites.   These 
field  measurements  were  obtained  using  inductance  coils  that  measure  ballast 
strain  and  optical  surveying  from  benchmarks. 

The  report  will  conclude  with  a  discussion  of  the  various  factors  that 
have  not  been  accounted  for  in  the  deformation  prediction  model,  and  the  im- 
plications of  these  factors  on  track  performance. 
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CHAPTER  3.   SITE  DESCRIPTIONS 


3.1  General  Background 

Four  revenue  service  track  locations  and  the  experimental  FAST  track  in 
Pueblo,  Colorado  were  selected  as  test  areas  for  this  research  project.   The 
revenue  service  sites  included  three  locations  that  contained  concrete  cross 
ties  and  a  control  section  having  wood  cross  ties.   The  wood  tie  control  sec- 
tion and  one  concrete  tie  section  are  located  at  Leeds  near  Streator,  Illi- 
nois, in  the  north-central  section  of  the  state.   The  remaining  two  concrete 
tie  test  sections  are  located  at  Aberdeen,  Maryland  and  Lorraine  near  Richmond, 
Virginia.   A  more  detailed  description  of  each  of  these  four  sites  was  given  in 
Refs.  [1,2,3,4,5].   The  FAST  track  contained  both  wood  and  concrete  tie  test 
sections. 

The  FAST  track  sections  17,  18  and  19  were  instrumented  in  1976,  when  the 
track  was  originally  constructed,  in  order  to  monitor  track  foundation  defor- 
mations and  stresses.   The  track  conditions  for  these  sections  are  given  in 
Table  la.   The  measurements  made  in  these  track  sections  are  used  in  portions 
of  the  correlation  study. 

Part  of  the  FAST  track  was  reconstructed  after  400  million  gross  tons 
(MGT)  traffic  had  accumulated  over  the  original  track  at  the  site.   The  gene- 
ral rebuild  procedure  for  these  test  sections  consisted  of: 

a)  Undercutting  the  existing  track. 

b)  Removal  of  the  existing  track  structure. 

c)  Subgrade  preparation,  including  grading  and  compaction  where  necessary. 

d)  Building  new  skeletonized  track  on  the  subgrade. 

e)  Dumping  and  tamping  ballast  and  raising  the  track. 

f)  Regulating  the  ballast,  followed  by  surfacing  and  lining  operations. 
Further  details  regarding  the  rebuild  procedures  can  be  found  in  Ref .  [6]. 

Data  for  certain  ballast  physical  state  tests  were  available  from  the  re- 
built sections  3  and  22.  The  track  conditions  for  these  sections  are  given  in 
Table  lb.   The  concrete  ties  used  in  section  22  were  the  RT7-SS2  type. 

3.2  Track  Structure 

The  track  structure  for  railroads  has  several  mechanical  components,  the 
most  obvious  of  which  are  the  cross  ties  and  rails.   Rails  are  generally  fas- 
tened to  the  wood  ties  using  the  conventional  tie  plate  and  cut  spike  system. 
Rails  are  attached  to  the  concrete  ties  using  a  spring-type  clip  fastener  sys- 
tem.  Tie  pads  support  the  rail  on  concrete  ties,  and  electrical  insulators 
separate  the  steel  rail  from  the  steel  rail  clips.   Table  2  lists  the  main 
track  structural  parameters  for  FAST  and  the  revenue  service  test  sections. 
Some  of  the  structural  characteristics  which  differ  among  the  test  sections 
are:   a)  tie  type  and  dimensions,  b)  tie  spacing,  c)  tie  weight,  d)  rail  sec- 
tion, e)  rail  fastener,  and  f)  pad  type. 

3.3  Ballast  Materials  and  Physical  States 

A  description  of  the  ballast  materials,  along  with  the  results  of  the 
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Table  1.  FAST  Track  Section  Conditions 


a)  Original  Track 


Tie 

Track 

Track 

Ballast 

Ballast 

Tie 

Spacing 

Section 

Geometry 

Type 

Depth  (in.) 

Type 

(in.) 

17C 

Curved 

Granite 

15 

Concrete 

24.0 

17E 

Tangent 

Granite 

15 

Concrete 

24.0 

18A 

Tangent 

Granite 

21 

Wood 

19.5 

18B 

Tangent 

Granite 

15 

Wood 

19.5 

20B 

Tangent 

Limestone 

15 

Wood 

19.5 

20G 

Tangent 

Trap rock 

15 

Wood 

19.5 

b)  Reconstructed  Track 


3 

Curved 

Granite* 

20 

Wood 

19.5 

17ABCDE 

Curved 

Granite 

16 

Concrete 

24.0 

22A 

Tangent 

Traprock 

15 

Wood 

19.5 

22B 

Tangent 

Traprock 

15 

Concrete 

24.0 

*Several  types  of  granite  were  used, 
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ballast  physical  state  tests  was  presented  in  Ref.  [1].  A  summary  of  the  re- 
sults and  conclusions  will  be  presented  here  where  relevant  to  the  prediction 
methodology. 

The  physical  state  tests  used  to  characterize  the  ballast  materials  were 
ballast  density,  plate  load, and  lateral  tie  push.   These  tests  were  done  be- 
fore and  after  the  track  was  raised  and  tamped  so  that  quantitative  assess- 
ments of  the  change  in  ballast  physical  state  could  be  made. 

The  ballast  density  tests  were  done  in  the  cribs  and  under  the  ties  as 
close  to  the  rail  seat  as  possible.   The  plate  load  tests  were  done  at  similar 
locations.   The  lateral  tie  push  tests  were  done  on  single  ties  located  away 
from  the  density  and  plate  load  test  areas. 

The  ballast  density  test  procedure  has  been  described  in  Refs.  [7  and  8]. 
The  test  results  showed  that  the  maintenance  operation  decreased  the  relative 
compaction  of  the  ballast  under  the  ties  by  10-12%  at  the  sites  where  a  large 
raise  was  applied.   Large  raises  of  1-1/2  to  2  in.  were  given  to  the  Lorraine 
and  Leeds  wood  tie  section.   The  pre-  and  post-maintenance  ballast  densities 
were  the  same  at  the  Aberdeen  site,  where  the  raise  was  only  about  0.1  in. 

Plate  load  tests  on  the  ballasts  were  done  before  and  after  maintenance 
at  selected  locations  in  the  cribs  and  under  the  ties.   Descriptions  of  the 
plate  load  test  procedures  and  apparatus  can  be  found  in  Refs.  [8,9,10],   The 
test  results  showed  that  the  largest  plate  bearing  resistances  were  under  the 
ties  near  the  rail  seat  areas.   The  pre-maintenance  values  were  greatest  at 
the  Lorraine  site,  possibly  because  of  the  high  degree  of  fouling  and  the 
greatest  accumulation  of  tonnage  since  previous  maintenance. 

The  post-maintenance  values  of  plate  bearing  index  were  very  similar  for 
all  of  the  sites,  with  the  exception  of  Aberdeen.   The  Aberdeen  site  showed 
the  smallest  change  in  plate  bearing  index  at  the  under-tie  location  as  a  re- 
sult of  the  small  raise.   The  largest  decrease  in  plate  bearing  index  due  to 
maintenance  was  also  at  the  Lorraine  site.   This  site  had  the  largest  raise, 
about  2  in. ,  which  would  significantly  disturb  the  ballast  structure,  particu- 
larly when  the  ballast  is  fouled.   This  was  consistent  with  the  observation 
that  the  Lorraine  site  also  showed  the  largest  decrease  in  relative  compaction 
at  the  under-tie  location.   For  all  of  the  sites  where  significant  trackbed 
disturbance  under  the  ties  resulted  from  raising  and  tamping,  the  plate  load 
tests  showed  no  significant  difference  in  ballast  stiffness,  independent  of 
ballast  type. 

Values  of  the  plate  bearing  index  were  also  measured  under  the  center  of 
the  tie.   Again,  the  Lorraine  site  had  the  highest  average  value  before  the 
maintenance.   The  Aberdeen  values  at  this  location  showed  no  change  resulting 
from  maintenance,  but  do  show  that  the  values  under  the  center  of  the  ties 
prior  to  maintenance  were  only  1/3  of  the  value  measured  under  the  tie  at  the 
rail  seat.   The  plate  load  tests  for  all  of  the  revenue  sites  showed  that  the 
ballast  stiffnesses  were  always  greater  under  the  tie  at  the  rail  seat  than 
under  the  tie  center  prior  to  maintenance.   After  the  maintenance  operations, 
the  ballast  stiffnesses  under  the  tie  were  about  the  same  at  both  locations, 
although  still  slightly  greater  near  the  rail  seat  areas. 

3.4  Ballast  Resilient  Behavior 

General  Concepts 

The  behavior  of  a  granular  material  under  repeated  loading  conditions  has 
generally  been  recognized  as  being  nonlinearly  stress-dependent,  and  very  dif- 
ferent from  the  behavior  of  the  same  soil  under  monotonic  loading  conditions. 
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This  type  of  behavior  is  shown  in  Fig.  1.   During  the  primary  loading  of  a 
granular  soil  such  as  ballast,  the  strain  develops  rapidly  and  is  only  par- 
tially recovered  upon  unloading.   Each  additional  cycle  contributes  another 
increment  of  plastic  or  permanent  strain.   The  magnitude  of  the  increment 
generally  decreases  with  the  number  of  cycles.   The  difference  between  the 
maximum  strain  under  peak  load  and  the  permanent  strain  after  unloading  for 
each  cycle  is  the  resilient  or  recoverable  strain.   The  resilient  modulus  of 
materials  is  defined  as  the  repeated  deviator  stress  divided  by  the  recovera- 
ble strain.   The  amount  of  resilient  recoverable  strain  generally  decreases 
with  number  of  cycles. 

Research  has  shown  that  under  moderate  levels  of  repeated  load,  the  re- 
silient strains  become  approximately  constant  after  some  number  of  cycles, 
and  the  material  behaves  elastically.   Kalcheff  and  Hicks  [11]  have  indicated 
that  only  a  few  hundred  cycles  are  necessary  for  this  stabilization  to  occur. 
However,  Brown  [12]  and  Morgan  [13]  have  shown  that  several  thousand  load 
repetitions  are  necessary  to  reach  constant  values  of  resilient  modulus. 

The  magnitude  of  the  resilient  modulus  is  very  much  stress-state  depen- 
dent.  Several  studies  [14,15,16]  have  shown  that  the  resilient  response  of 
unbound  granular  materials  greatly  increases  as  the  confining  pressure  in- 
creases and  is  affected  to  a  much  smaller  extent  by  the  magnitude  of  repeated 
deviator  stress.   This  has  led  to  the  conclusion  that  the  resilient  modulus 
can  be  related  to  the  bulk  stress  state  using  the  relationship 

_K2 

Er  =  KjB     ,  (1) 

where  E  =  resilient  modulus, 
r 

K..  ,K~  =  soil  constants  determined  from  laboratory  test  results,  and 

0  =  bulk  stress  or  the  first  invariant  of  stress. 
The  bulk  effective  stress  is  defined  by 


6  "  °1  +  °2  +  °3 

=   O     +  20     for  triaxial  compression  tests, 

where  O   ,  0  ,  O  =  major,  intermediate  and  minor  effective  principal  stresses, 
1   ^   -*   respectively. 

This  formulation  has  been  widely  used  [17,18,19]  to  characterize  the  resilient 
modulus  of  granular  materials. 

The  bulk  stress  formulation  (Eq.  1)  has  been  preferred  for  ballast  mater- 
ial by  Robnett,  et  al.  [20]  and  Knutson  [21].   Alva-Hurtado  [22]  has  investi- 
gated the  resilient  response  of  a  granite  railroad  ballast  and  determined 
that  an  arithmetic  relationship  was  equally  valid  when  relating  bulk  stress 
to  resilient  modulus.   The  relationship  developed  by  Alva-Hurtado  also 
showed  that  the  initial  compaction  state  affected  the  resilient  modulus. 

Although  it  has  been  shown  that  granular  materials  exhibit  high  values 
of  resilient  modulus  even  when  the  applied  deviator  stress  level  is  near  fail- 
ure [23],  the  use  of  bulk  stress  cannot  allow  any  distinction  to  be  made  be- 
tween the  effects  of  deviator  stress  or  cell  pressure.   Various  combinations 
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of  cell  pressure  and  deviator  stress  are  assumed  to  achieve  the  same  bulk, 
stress,  i.e.,  low  cell  pressure  -  high  deviator  stress  or  high  cell  pressure  - 
low  deviator  stress  or  medium  cell  pressure  -  medium  deviator  stress  combina- 
tions, all  can  result  in  the  same  bulk  stress. 

Staged  Testing 

In  order  to  study  the  resilient  behavior  of  railroad  ballast  under  more 
general  loading  conditions   that  would  include  partial  unloading  and  shear 
stress  reversal,  two  triaxial  test  series  were  performed.   The  first  series 
of  triaxial  tests  did  not  involve  shear  stress  reversal  so  the  major  principal 
stress  was  always  the  vertical  stress.   These  tests  were  called  staged  tests, 
since  a  sample  was  subjected  to  several  stress  paths  involving  full  and  par- 
tial unloadings  and  various  combinations  of  maximum  and  minimum  positive  shear 
stress.   Previous  research  by  Hicks  [14]  and  Allen  [24]  has  indicated  that  one 
specimen  can  be  used  to  measure  the  resilient  response  of  granular  materials 
over  a  range  of  stress  levels,  and  that  the  sequence  of  applied  stresses  did 
not  significantly  affect  the  measured  resilient  response. 

The  ballast  material  used  for  both  test  series  was  the  Wyoming  granite 
originally  used  at  FAST  in  sections  17  and  18.   This  material  was  chosen  for 
several  reasons.   The  static  compression  stress-strain  behavior  was  well  docu- 
mented [22,25],  data  and  bulk  stress  formulations  for  the  resilient  response 
had  been  generated  by  other  researchers  [17,22],  and  reproduceable  sample 
preparation  and  compaction  procedures  had  been  established. 

The  gradation  curve  for  the  granite  ballast  is  shown  in  Fig.  2.  This 
ballast  can  be  classified  as  an  AREA  Number  5  material.  The  specification 
ranges  for  the  AREA  Number  5  ballast  are  shown  in  Fig.  2. 

All  of  the  staged  cyclic  tests  were  done  on  compacted  samples  of  ballast. 
The  rubber-tipped  impact  hammer  compaction  method  developed  by  Yoo,  et  al. 
[7]  was  chosen  to  prepare  the  laboratory  triaxial  samples.   The  method  con- 
sisted of  placing  layers  of  ballast  into  a  specially  designed  split  compaction 
mold  that  was  clamped  together  and  placed  around  the  triaxial  cell  base  pedes- 
tal.  The  layer  was  then  subjected  to  twenty  blows  from  the  rubber- tipped  im- 
pact hammer.   This  compactive  effort  resulted  in  reproduceable  density  states. 

After  each  sample  was  compacted  in  the  mold,  the  sample  was  aligned  with 
the  triaxial  system  loading  piston  and  the  mold  was  removed.   A  small  vacuum 
was  applied  to  the  inside  of  the  sample  prior  to  removal  of  the  mold  in  order 
to  give  the  sample  enough  strength  so  that  it  would  not  collapse  during  fur- 
ther test  preparation.   The  ballast  was  tested  dry  to  provide  for  more  rapid 
sample  preparation  and  ease  of  testing.   The  effective  isotropic  consolida- 
tion pressures  were  achieved  by  increasing  the  vacuum  applied  to  the  inside  of 
the  sample,  eliminating  the  need  for  an  exterior  fluid  pressure  surrounding 
the  triaxial  sample  cylinder. 

The  specific  initial  test  conditions,  sample  dimensions,  and  dry  densi- 
ties for  the  staged  tests  are  given  in  Table  3,  along  with  the  results  from 
tests  involving  shear  stress  reversal.   As  can  be  seen  from  the  results 
given  in  Table  3,  the  compaction  procedures  gave  reproduceable  densities.   The 
average  dry  density  was  100.3  pounds  per  cubic  foot  (pcf)  and  had  a  coefficient  of 
variation  of  only  0.7%.   The  average  initial  void  ratio,  eo,  was  0.66. 

As  previously  stated,  the  purpose  of  the  staged  cyclic  tests  was  to 
investigate  the  relationship  between  the  applied  stress  level  and  resilient 
modulus  for  both  full  unloading  and  partial  unloading.   Full  unloading  for 
these  tests  will  result  from  complete  removal  of  any  applied  deviator  stress, 
O  -O    .      Partial  unloading  will  result  from  cycling  between  two  levels  of  devia- 
tor stress  without  allowing  complete  unloading. 
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Confining  pressures  of  3,  6  and  12  psi  were  selected  for  these  staged 
tests.   For  each  confining  pressure,  four  levels  of  applied  shear  stress  (q) 
were  used.   The  various  combinations  of  cell  pressure,  maximum  shear  stress, 
and  minimum  shear  stress  used  for  the  test  stress  paths  are  given  in  Table  4. 
Table  4  also  gives  the  maximum  bulk  stress,  G,  for  each  test,  along  with  the 
stress  ratio,  q    /qf>  where  q   is  the  static  compressive  strength  of  the  bal- 
last  for  a  given  effective  cell  pressure. 

Four  individual  staged  tests  were  done  at  each  confining  stress.   In  each 
test,  the  order  of  stress  application  was  varied  to  investigate  any  changes  in 
resilient  modulus  that  might  result  from  variations  in  stress  sequence.   Per- 
manent deformation  characteristics  for  each  stress  level  and  sequence  were 
also  analyzed.   These  results  will  be  discussed  in  another  section. 

The  order  of  applied  stress  paths  for  each  of  the  tests  is  given  in 
Table  5.   Each  stress  path  was  repeated  for  1000  cycles  of  loading  and  unload- 
ing and  the  recoverable  strain  between  the  maximum  and  minimum  applied  shear 
stress  was  measured.   The  recoverable  strains  during  cycling  along  any  of  the 
stress  paths  tended  to  decrease  as  the  number  of  cycles  increased,  but  had 
generally  reached  an  equilibrium  value  after  about  500  cycles  of  the  stress 
path  had  been  applied.   The  resilient  modulus,  E  ,  resulting  from  any  of  the 
stress  paths  was  calculated  by  dividing  the        resilient  strain  into  twice 
the  cyclic  shear  stress,  Aq,  according  to 

e   =^a     ,  (2) 

re 
r 

where  Aq  =  difference  between  maximum  and  minimum  applied  shear  stress,  and 

£   =  resilient  strain, 
r 

The  averaged  results  of  the  resilient  moduli  for  each  of  the  individual 
stress  paths  are  given  in  Table  6,  along  with  the  corresponding  average  resi- 
lient strains.   The  resilient  moduli  values  given  in  Table  6  are  plotted 
against  the  maximum  bulk  stress  for  each  stress  path  in  Fig.  3.   There  is  no 
distinct  relationship  between  bulk  stress  and  resilient  modulus  evident  when 
moduli  values  measured  from  partial  unloading  stress  paths  are  included.   A 
check  was  also  made  to  see  if  the  data  fit  the  common  log-log  f ormulation. 
Again,  there  was  good  correlation  between  log  of  E  and  log  of  9  if  only  the 
results  from  fully  unloaded  stress  paths  were  used. 

Figure  4  shows  the  comparison  between  resilient  moduli  values  obtained 
by  Alva-Hurtado  [22]  and  those  obtained  in  this  study.   The  arithmetic  regres- 
sion equations  for  these  two  sets  of  data  are  given  in  Table  7a.   Although 
both  sets  of  moduli  measurements  were  made  on  the  same  granite  material  at  the 
same  density  state,  the  results  obtained  by  Alva-Hurtado  are  higher  than  those 
obtained  in  this  study.   The  reason  for  this  is  not  known,  but  it  could  be 
from  differences  in  detailed  test  conditions  or  from  changes  in  ballast  grada- 
tion and  angularity  caused  by  repeated  use  of  the  same  material.   However,  the 
slopes  of  the  two  E  vs  0  curves  shown  in  Fig.  4  are  the  same.   The  coeffi- 
cients for  the       log-log  regression  equations  for  the  resilient  moduli- 
bulk  stress  formulation  from  these  sets  of  data  are  also  given  in  Table  7b  for 
comparison. 

Shear  Stress  Reversal  Testing 

The  resilient  behavior  of  ballast  under  one-way  repeated  loading  in  which 
the  full  deviator  stress  was  removed  was  shown  to  be  related  to  the  maximum 
bulk  stress  on  the  material.   Partial  unloading  resilient  moduli  are  not 
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Table  4.   Stress  Levels  for  Staged  Tests 


Stress          3c 

o 

1 

q 

Aq 

e 

Max 

nnax 

Min 

Max 

Min 

Max 

Path        (psi) 

(psi) 
3.00 

(psi) 
6.12 

(psi) 
0.00 

(psi) 
1.56 

(psi) 
1.56 

(psi) 
12.12 

qf 

OA            3.00 

0.150 

OB 

3.00 

9.26 

0.00 

3.13 

3.13 

15.26 

0.301 

OC 

3.00 

12.38 

0.00 

4.69 

4.69 

18.38 

0.451 

OD 

3.00 

15.51 

0.00 

6.26 

6.26 

21.51 

0.602 

AB 

6.12 

9.26 

1.56 

3.13 

1.56 

15.26 

0.301 

AC 

6.12 

12.38 

1.56 

4.69 

3.13 

18.38 

0.451 

AD 

6.12 

15.51 

1.56 

6.26 

4.69 

21.51 

0.602 

BC 

9.26 

12.38 

3.13 

4.69 

1.56 

18.38 

0.451 

BD 

9.26 

15.51 

3.13 

6.26 

3.13 

21.51 

0.602 

CD 

12.38 

15.51 

4.69 

6.26 

1.56 

21.51 

0.602 

O'a'          6.00 

6.00 

12.00 

0.00 

3.00 

3.00 

24.00 

0.150 

0  B 

6.00 

18.00 

0.00 

6.00 

6.00 

30.00 

0.301 

O'c' 

6.00 

24.00 

0.00 

9.00 

9.00 

36.00 

0.451 

O'D' 

6.00 

30.00 

0.00 

12.00 

12.00 

42.00 

0.602 

A"B' 

12.00 

18.00 

3.00 

6.00 

3.00 

30.00 

0.301 

A'c' 

12.00 

24.00 

3.00 

9.00 

6.00 

36.00 

0.451 

A  D 

12.00 

30.00 

3.00 

12.00 

9.00 

42.00 

0.602 

B  C 

18.00 

24.00 

6.00 

9.00 

3.00 

36.00 

0.451 

B  D 

18.00 

30.00 

6.00 

12.00 

6.00 

42.00 

0.602 

C'D' 

24.00 

30.00 

9.00 

12.00 

3.00 

42.00 

0.602 

0"A"   12.00 

12.00 

18.00 

0.00 

3.00 

3.00 

42.00 

0.084 

0~B~ 

12.00 

24.00 

0.00 

6.00 

6.00 

48.00 

0.168 

0~C~ 

12.00 

30.00 

0.00 

9.00 

9.00 

54.00 

0.252 

0~D" 

12.00 

36.00 

0.00 

12.00 

12.00 

60.00 

0.335 

A^B" 

18.00 

24.00 

3.00 

6.00 

3.00 

48.00 

0.168 

k"or* 

18.00 

30.00 

3.00 

9.00 

6.00 

54.00 

0.252 

a"d" 

18.00 

36.00 

3.00 

12.00 

9.00 

60.00 

0.335 

z"c" 

24.00 

30.00 

6.00 

9.00 

3.00 

54.00 

0.252 

B~D~ 

24.00 

36.00 

6.00 

12.00 

6.00 

60.00 

0.335 

C"D" 

30.00 

36.00 

9.00 

12.00 

3.00 

60.00 

0.335 
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Table  6.    Average  Resilient  Moduli  and  Strains  for  Staged  Tests 


2Aq/E 

_ 

r 

Q3c 

Stress 

Aq 

Avg.    Ef 

=  Avg.    £ 
r 

(psi) 

Path 
OA 

(psi) 
1.56 

(ksi) 
14.4 

(ye) 

3.0 

217 

OB 

3.13 

15.7 

399 

OC 

4.69 

16.8 

559 

OD 

6.26 

17.7 

709 

AB 

1.56 

25.4 

123 

AC 

3.13 

24.3 

252 

AD 

4.69 

24.6 

381 

BC 

1.56 

33.1 

94 

BD 

3.13 

30.5 

206 

CD 

1.56 

40.8 

76 

6.0 

0"A' 

3.00 

23.2 

259 

O'B' 

6.00 

24.0 

500 

0"C 

9.00 

25.6 

704 

O'D' 

12.00 

27.0 

888 

A'B' 

3.00 

38.8 

154 

A'C 

6.00 

37.8 

318 

A'D' 

9.00 

37.9 

475 

B'C* 

3.00 

53.3 

113 

B'D' 

6.00 

45.1 

266 

C'D" 

3.00 

57.5 

104 

12.0 

0"A" 

3.00 

33.0 

182 

0~B~ 

6.00 

32.8 

366 

o"c" 

9.00 

33.3 

540 

0"D~ 

12.00 

33.9 

707 

A~B" 

3.00 

50.3 

119 

k"Z" 

6.00 

43.6 

275 

A~D~ 

9.00 

44.0 

409 

B"C" 

3.00 

68.1 

88 

B~D~ 

6.00 

58.0 

207 

C"V" 

3.00 

84.9 

71 
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Table  7.    Regression  Equations  for  E  as  a  Function  of  9 


a)  Model:  E  =  K  +  K  0 


K, 


Staged  Tests 


9690 


449 


0.92 


Alva-Hurtado  [22] 
(N  =  10,000) 


22685 


425 


0.94 


b)  Model:  E  =  Ko0 
r    3 


K. 


Staged  Tests 


3160 


0.593 


0.95 


Alva-Hurtado  [22 
(N  -  10,000) 


5477 


0.535 


0.91 
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related  to  the  bulk  stress  alone,  and  thus  must  be  influenced  by  a  combination 
of  stress  parameters.   When  triaxial  samples  are  subjected  to  both  positive 
and  negative  cyclic  shear  stresses,  the  resilient  behavior  is  even  more  com- 
plicated.  Cycling  between  positive  and  negative  shear  stresses  causes  a  con- 
tinual reversal  of  the  principal  stress  directions  in  a  triaxial  sample.   For 
triaxial  compression,  the  final  stress  state  is  such  that  the  major  principal 
stress  acts  in  the  vertical  direction  on  a  horizontal  plane.   For  triaxial 
extension,  the  final  major  principal  stress  acts  in  the  horizontal  direction 
on  a  vertical  plane. 

A  "positive"  shear  stress,  q  =  (a  -a  )/2,  will  be  defined  as  the  shear 
stress  resulting  from  a  major  principal  stress  acting  within 

±45°  from  the  vertical.   Likewise,  a  "negative"  shear  stress  occurs  when  the 
major  principal  stress  is  acting  within  ±45°  from  the  horizontal.   Thus, 
these  triaxial  compressional  stresses  produce  positive  shear  stresses  and 
these  extensional  stresses  produce  negative  shear  stresses.   When  the  shear 
stresses  are  cycled  from  positive  to  negative,  the  orientation  of  the  major 
principal  stress  and  plane  is  abruptly  rotated  by  90°  as  the  stresses  pass 
through  the  isotropic  condition. 

The  sample  initial  dimensions,  densities  and  void  ratios  for  the  stress 
reversal  tests  have  been  given  in  Table  3.   The  test  conditions  for  these 
tests  are  given  in  Table  8.   For  these  tests,  O     will  mean  the  vertical  stress 
acting  on  the  sample,  which  can  be  the  major      principal  stress,  a  ,  or  the 
minor  principal  stress,  o   depending  on  the  overall  stress  condition.   Simi- 
larly, a  will  refer  to    the  horizontal  stress  acting  on  the  sample.   All 
test  samples  were  isotropically_  consolidated  prior  to  shear  application.   For 

the  isotropic  condition,  o  =o=o=g=g,   where  O   is  the  isotropic 

v    h    1    3    c         c 
consolidation  pressure.  '  Also  given  in  Table 

8  are  the  stress  ratios,  q/qf,  for  both  the  compressional_and  extensional  por- 
tions of  the  applied  stress  paths,  and  the  mean  stress,  p,  applied  to  the 
sample. 

For  triaxial  extension  tests,  the  static  shear  strength  is  defined  by 


c  cos4>  -  °,  f  sin$ 

" •  (3) 

1  +  sin4> 


lf 


where  o  f  is  the  major  principal  stress  at  failure.   For  these  triaxial  exten- 
sion     tests,  0=o=o    ,  since  the  cell  pressure  was  maintained  con- 

,     ,      It    hi    c  , 
stant  throughout  the  test. 

A  cycle  of  loading  for  the  stress  reversal  tests  consisted  of  first  ap- 
plying the  axial  compressive  stress  followed  by  unloading  to  the  extensional 
stress,  and  then  reloading  to  isotropic  conditions.   Loading  is  defined  as 
any  change  of  stress  state  where  the  mean  stress  increases,  that  is,  p  = 
(a  +a  )/2  increases.   Unloading  is  defined  as  a  decrease  in  p.   During  testing, 

the  resilient  strains  resulting  from  loading  to  q    and  unloading  to 

nmax 


\ 


in 


were  monitored.   The  resilient  moduli  calculations     were  based  on 

=  _Vn  (4) 

re     e 

re 

E   =i^ax  (5) 


re     £ 
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2(q    -q  . 
_    max  nmin 

Ert  "  E   -  £       '  (6) 

re    re 

where  E   ,  E   and  E  ,  are  the  extension,  compression  and  total  resilient 
re   re      rt    ,  1  . 
moduli, 

£   and  e   are  the  resilient  strains  resulting  from  extension  and 

compression,  and 

q    and  q  .   are  the  maximum  and  minimum  shear  stresses, 
max      min 

The  measured  resilient  moduli  for  these  tests  are  listed  in  Table  9  for  var- 
ious cycle  numbers.   These  tabulated  values  show  that  the  extensional  resi- 
lient modulus  is  always  less  than  the  compressional  modulus,  which  results  in 
the  total  resilient  modulus  having  a  value  intermediate  to  the  compressional 
and  extensional  moduli.   The  general  trend  for  these  tests  is  that  the  moduli 
increase  as  number  of  cycles  increases.   In  contrast,  the  results  from  the 
staged  tests  showed  that  the  moduli  reached  approximately  constant  values 
after  a  fewer  number  of  cycles. 

Resilient  Modulus  Formulation 

After  many  load-unload  cycles,  the  behavior  of  railroad  ballast  becomes 
relatively  elastic,  in  that  the  resilient  strains  become  approximately  con- 
stant.  The  elastic  behavior  is  nonlinear,  however. 

The  resilient  strains  from  the  loading  portions  of  the  staged  tests  and 
the  peak  points  of  the  shear  stress  reversal  tests  are  given  in  Table  10. 
These  points  will  be  used  to  determine  the  resilient  modulus  formulation  for 
situations  involving  full  unloading  and  with  shear  stress  reversal.   The  re- 
sults from  the  stress  reversal  test  at  9  psi  effective  confining  pressure  are 
not  included  because  equipment  malfunction  during  testing  invalidated  some  of 
the  test  results. 

Figure  5  shows  the  relationships  between  the  applied  shear  stress  levels 
and  the  vertical  resilient  strains.   The  resilient  strains  increase  as  cell 
pressure  decreases,  for  a  given  shear  stress  level,  for  both  extension  and 
compression. 

The  compressional  resilient  strains  in  Fig.  5  for  each  cell  pressure 
were  fit  to  an  equation  of  the  form 

q  =  A(U£r)B    ,  (7) 

using  the  data  listed  in  Table  10.   Since  intermediate  data  were  not  availa- 
ble for  the  extensional  resilient  strain  curves,  straight  lines  from  the  ori- 
gin through  the  peak  extensional  resilient  strains  were  assumed.   The  resi- 
lient extensional  strains  were  related  to  the  negative  shear  stress  by 

q  =  a  (ye)   .  (8) 

The  coefficients  A,  B  and  a  for  each  cell  pressure  were  subsequently  fit  to  a 
log-log  relationship  as  a  function  of  cell  pressure.   Good  agreement  between 
the  observed  data  and  the  assumed  forms  of  the  regression  equations  was  found 
for  the  fitted  data.   The  solid  lines  shown  on  Fig.  5  were  based  on  the  regres- 
sion equations,  and  are  in  close  agreement  with  the  observations. 

Table  11  lists  the  assumed  form  of  each  regression  equation  and  the  de- 
rived coefficients  and  coefficients  of  determination  for  the  regressions.   The 
resilient  strain  relationships  given  in  Table  11  were  used  to  define  contours 
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Table  10.       Resilient   Strains    for  Staged  Tests   and  Stress 
Reversal  Tests 

a)    Staged  Tests 

<L  e 

3                                                                 q  r 

(psi)                                                            (psi)  (lie) 

3.0                                                           1.56  217 

3.13  399 

4.69  560 

6.26  709 

6.0                                                           3.00  259 

6.00  500 

9.00  704 

12.00  888 

12.0                                                           3.00  182 

6.00  366 

9.00  540 

12.00  707 

b)    Stress  Reversal  Tests 

a0  e 

3                                                                 q  r 

(psi)                                                         (psi)  (ye) 

6.0                                                        -1.83  -467 

2.99  242 

12.0                                                        -3.01  -318 

-2.99  -264 

-2.99  -293 

2.99  182 

5.98  320 
9.11  477 

15.0                                                        -3.01  -184 

8.99  415 
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Table  11.   Regression  Results  for  Compression  and 
Extension  Resilient  Strains 

a)  Compression 

■Q 

Form  of  Equation:   q  ■  A(ye  ) 


w3 

(psi) 

A 

B 

3.0 

2.486  x  10"3 

1.171 

6.0 

7.680  x  10"3 

1.079 

12.0 

13.914  x  10"3 

1.036 

1.00 
1.00 
0.99 


Equations  for  Parameters  A,  B: 

A  =  c(a3)D 


C  =  8.647  x  10 


-4 


B  =  E(a3)F 


E  =  1.2816 


D  =  1.145 


-0.0883 


r  =  0.98 


r  =  0.96 


b)  Extension 


Form  of  Equation: 


q  =  a(ye  ) 
r 


u3 
(psi) 

6.0 

12.0 

15.0 


a 


3.921  x  10 


-3 


1.023  x  10 


-3 


1.639  x  10 


-3 


Equation  for  a 


-  nY 


a  =  S(a3) 


3  =  2.535  x  10 
'  y  =  1.518 
r2  =  0.99 


-4 
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of  equal  resilient  strain  as  related  to  p  =  0  +  0,  2   and  q  =  0  -0 ,/2  for  full 

-   , .  ml       ....        .      V   h  vh 

unloading  cases.   These  resilient  strain  contours 

are  shown  in  Fig.  6.   This  figure  shows  that  for  a  given  effective  confining 
pressure,  a  positive  shear  stress  results  in  smaller  resilient  strains  than 
the  application  of  the  same  magnitude  negative  shear  stress.   This  means  that 
the  compressional  resilient  modulus  will  be  larger  than  the  extensional  resi- 
lient modulus  for  isotropic  tests  cycled  between  ±  q. 

The  procedure  for  calculating  the  resilient  modulus,  E^ ,  for  compres- 
sion loading  from  q  =  0  to  a  specified  value  of  q 
loading  to  q  =  0  is  as  follows: 

1)  The  expression  relating  shear  stress,  q 

as  a  function  of  effective  confining  pressure, 

to 

r 

-F 


re 


max 


max  — 


followed  by  full  un- 
to resilient  strain,  £ 


'3' 


can  be  rearranged 


re 


U£ 


rc 


-D 


max  3 


(9) 


where  C,  D,  E  and  F  are  the  coefficients  listed  in  Table  11.   Substitu- 
tion of  the  coefficients  results  in 

-  0.088] 


ye 


rc 


-1.145 


q     °-5 

max  3 


,-4 


1.282 


(10) 


(8.647  x  10 
2)  The  compressional  resilient  modulus,  E 

2  v 


rc 


is  then 


lax 


rc 


ue   x  10 
rc 


-6 


(11) 


The  procedure  for  calculation  of  the  extensional  resilient  modulus,  E   , 

re 
for  repeated  loadings  from  q  =  0  to  q.   is  as  follows: 

1)  The  expressions  relating  shear  stress,  q  .  ,  to  resilient  extensi_onal 
strains,  £  ,  as  a  function  of  effective       confining  pressure,  cr  , 
can  be  rearranged  to 


-Y 


ye. 


re 


Tnin  3 

3 


where  3  and  Y  are  the  coefficients  listed  in  Table  11. 
the  coefficients  results  in 


(12) 
Substitution  of 


-1.518 


ye 


mm  3 
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(13) 


2.535  x  10 
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2  q  . 
E   = 5in 

re   U£    x  10"6 
re 

For  tests  where  the  shear  stress  is  reversed,  the  repeated  stresses  are 

cycled  between  q    and  q_ .  .   The  total  resilient  modulus,  E   ,  would  be  cal- 

J ,  nmax     ^nin  ,  , .  rt 

culated  as  follows: 

1)  Calculate  ye    for  q    as  given  in  Eq.  9. 

re     nmax  n 

2)  Calculate  ye    for  q  .   as  given  in  Eq.  13. 

re     ^nin    °         n 

3)  Calculate  E  ,  by 

rt 

2(q    -q  .  ) 

max  mm  ,„  „. 

Ert  =  =6    '  (15) 

ru   (ye  -ye   )  x  10 
re   re 

Resilient  Modulus  Validation 

The  formulation  presented  for  full  unloading,  with  and  without  shear 
stress  reversal,  will  be  used  for  characterization  of  the  resilient  behavior 
of  the  ballast  materials. 

This  method  for  calculation  of  resilient  moduli  was  used  to  compute  the 
full  unloading  moduli  from  the  staged  and  stress  reversal  tests  on  the  compac- 
ted granite  ballast.   The  predicted  and  measured  resilient  moduli  values  for 
these  tests  are  given  in  Table  12. 

Repeated  load  tests  were  performed  on  compacted  samples  of  ballast  recov- 
ered from  the  Leeds,  Illinois, wood  and  concrete  tie  sections  and  the  Lorraine, 
Virginia,  and  Aberdeen,  Maryland,  concrete  tie  sections.   These  ballasts  and 
their  properties  were  previously  described  in  Volume  1  [1].   Details  regarding 
the  test  procedure  can  be  found  in  Ref.  [26].   The  measured  and  predicted  re- 
silient moduli  for  the  field  site  ballasts  after  about  10,000  load  cycles  are 
given  in  Table  13. 

The  comparisons  between  all  the  predicted  and  measured  resilient  moduli 
given  in  Tables  12  and  13  are  shown  in  Fig.  7.   Agreement  is  good  for  the 
staged  and  shear  stress  reversal  tests.   The  limited  data  on  resilient  moduli 
from  the  field  sites  show  that  the  field  site  ballasts  have  resilient  moduli 
similar  to  the  granite  ballast  moduli  measured  .in  the  staged  and  stress  rever- 
sal tests.   This  suggests  that  ballast  type  is  not  an  important  consideration 
with  respect  to  resilient  modulus. 

3.5  Ballast  Inelastic  Behavior 

Introduction 

Repeated  load  triaxial  tests  were  done  on  compacted  ballast  samples  to 
investigate  the  amount  of  permanent  deformation  that  accumulated  as  a  func- 
tion of  stress  level  and  number  of  constant-amplitude  load  applications. 
Variable-amplitude  repeated  load  tests  with  both  full  and  partial  unloading 
(staged  tests) ,  and  constant-amplitude  repeated  load  tests  with  shear  stress 
reversal  were  also  done  on  the  FAST  granite  ballast. 

The  stress  conditions  for  the  constant-amplitude  repeated  load  tests  have 
been  given  in  Table  13.   The  stress  conditions  and  load  sequences  for  the 
staged  tests  on  the  granite  ballast  have  been  given  in  Tables  4  and  5.   The 
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Table  12.   Resilient  Moduli  for  Granite  Ballast 
Staged  and  Stress  Reversal  Comparisons 


a3 

%n±n 
(psi) 

nnax 
(psi) 

Resilient  Modi 

llus,   E      (ksi) 

(psi) 

Measured 

Predicted 

3.0 

0.0 

1.56 

14.4 

14.6 

3.0 

0.0 

3.13 

15.7 

16.2 

3.0 

0.0 

4.69 

16.8 

17.1 

3.0 

0.0 

6.26 

16.7 

17.8 

6.0 

0.0 

3.00 

23.2 

22.8 

6.0 

0.0 

6.00 

24.0 

24.3 

6.0 

0.0 

9.00 

25.6 

25.1 

6.0 

0.0 

12.00 

27.0 

25.8 

12.0 

0.0 

3.00 

33.0 

34.8 

12.0 

0.0 

6.00 

32.8 

35.5 

12.0 

0.0 

9.00 

33.0 

35.9 

12.0 

0.0 

12.00 

33.9 

36.2 

6.0 

-1.8 

3.00 

13.6 

13.1 

9.0 

-3.0 

9.00 

19.8 

23.8 

12.0 

-3.0 

3.00 

24.0 

27.0 

12.0 

-3.0 

6.00 

30.7 

29.5 

12.0 

-3.0 

9.00 

31.4 

31.0 

15.0 

-3.0 

9.00 

40.1 

37.8 
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Table  13.   Resilient  Moduli  for  Field  Site  Ballasts 


Resilient  Modulus 


Site 


Leeds,  Wood 
(slag) 

Leeds,  Concrete 

(granite) 


Lorraine 

(limestone) 


Aberdeen 
(traprock) 


CT3 

qmin 
(psi) 

q 

max 

(psi) 

E 

r 

(ksi) 

(psi) 

Measured 

Predicted 

3.0 

0.0 

7.50 

23.0 

18.3 

12.0 

0.0 

7.50 

42.8 

35.7 

3.0 

0.0 

3.55 

17.7 

16.4 

3.0 

0.0 

7.50 

20.8 

18.3 

3.0 

0.0 

14.15 

17.0 

20.0 

12.0 

0.0 

3.75 

53.6 

35.0 

12.0 

0.0 

7.50 

38.8 

35.7 

12.0 

0.0 

15.00 

38.8 

36.4 

3.0 

0.0 

5.70 

27.6 

17.6 

3.0 

0.0 

7.50 

19.1 

18.3 

7.0 

0.0 

7.50 

27.8 

26.7 

12.0 

0.0 

7.50 

40.9 

35.7 

12.0 

0.0 

15.00 

43.1 

36.4 

3.0 

0.0 

3.00 

28.6 

16.1 

3.0 

0.0 

6.00 

33.8 

17.7 

12.0 

0.0 

3.75 

36.3 

35.0 

12.0 

0.0 

7.50 

41.3 

35.7 

12.0 

0.0 

15.00 

34.6 

36.4 
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test  conditions  for  the  shear  stress  reversal  tests  on  the  granite  ballast 
were  given  in  Table  8. 

In  addition  to  the  triaxial  tests,  experiments  were  conducted  using  a 
special  box  test  intended  to  simulate  the  field  loading  conditions.   The  de- 
formation properties  of  ballast  due  to  a  repeated  application  of  vertical 
loads  were  measured  in  this  device.   Lateral  stresses  were  also  measured  un- 
der both  loaded  and  unloaded  conditions,  and  residual  stress  build-up  in  the 
ballast  was  determined. 

Constant-Amplitude  Tests 

The  permanent  strains  resulting  from  the  repeated  loadings  of  the  reve- 
nue field  site  ballasts  have  been  previously  reported  by  Siller  [26].   The 
permanent  strains  developed  during  the  first  thousand  cycles  of  each  of  the 
staged  tests  are  given  in  Table  14,  and  will  be  considered  in  this  section. 
However,  the  staged  tests  were  carried  out  to  14000  cycles,  and  the  loading 
conditions  changed  every  1000  cycles.   The  permanent  strains  for  all  but  the 
first  stage  will  be  presented  later. 

Previous  research  on  the  permanent  strain  behavior  of  ballast  under  re- 
peated loading  conditions  has  shown  that  under  moderate  stress  levels,  a 
linear  relationship  exists  between  the  permanent  strain  at  the  Ntn  cycle  and 
the  log  of  the  number  of  cycles  [22,27,28,29].   This  relationship  has  been  ex- 
pressed by 

eN  =  E1(l  +  C  log  N)    ,  (16) 

where  E„  is  the  permanent  axial  strain  at  the  N   cycle, 

N 

Z      is  the  permanent  axial  strain  developed  due  to  the  first  cycle  of 
load,  N  =  1,  which  is  also  the  intercept  of  the  £  vs.  log  N  rela- 
tionship, 
C  is  a  soil  constant,  and 
N  is  the  number  of  applied  load  cycles. 
Past  research  has  indicated  that  C  -  0.2  for  all  ballasts  [22,27,28],  indepen- 
dent of  the  applied  stress  levels  or  compaction  state. 

The  ballast  constant,  C,  can  be  calculated  in  several  ways,  each  of  which 
should  give  the  same  result  if  the  relationship  given  in  Eq.  16  holds.   Four 
methods  for  calculating  C  are  as  follows: 

1.  Method  1  is  the  Office  of  Research  and  Experiments  (ORE)  method  [28], 
which  has  been  frequently  used  in  the  past  [22,27].  Details  of  the  ORE  method 
will  not  be  presented,  since  they  have  been  clearly  outlined  in  past  works. 

2.  Method  2  is  a  regression  of 

£N 
( 1)  vs.  log  N  to  calculate  C. 

€l 

3.  Method  3  determines  C  by  calculating 

£. 


1 
log  N 


for  all  N. 


4.   Method  4  uses  the  slope  of  the  semi-log  £   vs.  log  N  relationship  di- 
vided by  the  strain  at  the  first  cycle  to  calculate  C. 
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The  resulting  values  for  the  coefficient  C  for  each  of  the  ballasts  test- 
ed are  given  in  Table  15.   The  results  given  in  Table  15  show  that  all  of  the 
methods  used  were  in  very  close  agreement  for  the  revenue  site  ballasts  and 
for  the  staged  tests  done  on  the  FAST  granite. 

Since  the  four  methods  for  determination  of  the  coefficient  C  gave  consis- 
tent results  for  the  revenue  site  ballasts  and  the  granite  ballast  tested  most 
recently,  the  average  values  given  in  Table  15  were  taken  to  be  representative 
of  those  materials.   These  soil  constants  and  the  measured  first  cycle  strains 
were  used  to  calculate  strains  from  Eq.  16  for  comparison  with  the  measured 
strains  after  10,000  cycles  for  the  revenue  site  ballasts  and  after  1000  cyc- 
les for  the  FAST  granite.   As  indicated  in  Fig.  8,  the  agreement  is  very  good. 

Variable-Amplitude  Tests 

The  results  for  the  staged  tests  at  effective  confining  pressures  of  a   =  3, 
6  and  12  psi  are  shown  in  Figs.  9,  10  and  11,  respectively.   The  magnitude 
and  sequence  of  the  applied  stresses  used  in  these  variable-amplitude  tests 
are  given  in  Tables  4  and  5.   Several  very  clear  trends  were  evident  for  the 
staged  test  results.   They  were: 

1.  For  a  given  effective  cell  pressure,  the  permanent  strain  at  the 
first  cycle  increased  as  the  applied  vertical  deviator  stress  increased. 

2.  For  a  given  applied  vertical  deviator  stress,  the  permanent  strain  at 
the  first  cycle  decreased  as  the  effective  confining  pressure  increased. 

3.  When  the  deviator  stresses  were  increased  to  values  greater  than  any 
past  maximum  values,  the  permanent  strains  continued  to  increase. 

4.  When  the  deviator  stresses  were  reduced  to  values  less  than  any  past 
maximum  values,  negligible  changes  in  permanent  strain  resulted  with  addition- 
al cycling. 

5.  Partial  unloadings  resulted  in  negligible  increases  in  permanent 
strains  with  cycles. 

6.  For  any  given  effective  confining  pressure,  the  sequence  of  applied 
stresses  did  not  affect  the  final  values  of  permanent  strain  so  long  as  the 
total  number  of  cycles  at  each  stress  level  was  held  constant. 

Items  1  through  3  listed  above  are  well  known  and  easily  understood. 
Items  4  through  6  deserve  attention.   The  experimental  results  showed  that  the 
maximum  stress  history  of  a  sample  controls  future  stress-strain  response. 
These  ballast  test  results  imply  that  the  maximum  wheel  loads  for  a  given 
maintenance  life  cycle  will  control  the  amount  of  strain  accumulation. 

The  observation  that  maximum  load  controls  future  cyclic  deformations 
has  been  reported  for  other  ballast  materials  [21,27].   However,  past  re- 
search on  ballast  behavior  under  repeated  load  [27]  and  other  granular  mater- 
ials [30]  has  indicated  that  the  sequence  of  applied  loads  does  have  an  ef- 
fect on  the  total  permanent  strain  developed,  and  that  a  sequence  beginning 
with  high  loads  will  result  in  larger  final  permanent  strains.   This  was  not 
evident  from  the  staged  test  results  presented  here,  because  the  terminal 
permanent  strains  were  not  related  to  the  sequence  of  applied  stresses.   Con- 
sidering the  accepted  variability  of  sample  preparation  techniques  and  general 
ballast  behavior,  the  strains  measured  at  the  end  of  all  staged  test  series 
for  similar  effective  confining  stresses  were  in  very  close  agreement. 

The  small  amount  of  permanent  strain  accumulation  associated  with  partial 
unloadings  was  also  observed  by  Shenton  [27]  for  ballast  material. 

Shear  Stress  Reversal  Tests 

The  linear  relationship  between  permanent  vertical  strain  and  log  number 
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Figure  8.   Measured  versus  Calculated  Permanent  Strains  from  Triaxial 
Tests 
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of  cycles  was  evident  for  these  tests.   The  test  samples  all  continued  to  de- 
velop increasing  amounts  of  compressive  axial  strain  after  the  first  cycle. 
The  general  trends  noted  for  the  constant  repeated  stress  tests  with  full  un- 
loading were  also  observed  for  the  stress  reversal  tests.   As  effective  con- 
fining pressure  increased,  the  accumulated  permanent  strains  decreased  for 
constant  values  of  q    and  Aq.   For  constant  effective  confining  pressures 
and  equivalent  levels  of  negative  shear  stress,  the  amount  and  rate  of  accu- 
mulation of  permanent  strain  increased  as  the  compressive  shear  stress  levels 
increased. 

The  permanent  strain  developed  after  the  first  cycle  of  loading  for  tests 
involving  shear  stress  reversal  was  not  as  clearly  defined  as  for  the  one-way 
loading  tests.   The  strain  at  N  =  1  was  in  fact  extensional  (negative)  for  one 
of  the  stress  reversal  tests.   The  sample  did,  however,  still  tend  to  compress 
for  the  subsequent  cycles.   The  explanation  for  the  negative  axial  permanent 
strains  that  could  be  developed  for  these  stress  reversal  tests  is  clearly 
shown  in  Fig.  12. 

The  permanent  vertical  strain  developed  due  to  compression  loading  fol- 
lowed by  unloading  with  shear  stress  reversal  is  dependent  upon  not  only  the 
minimum  negative  shear  stress,  q.  ,  but  also  the  maximum  positive  shear 
stress.   For  the  same  level  of  q  .  ,  the  vertical  permanent  strain  at  the  end 
of  the  first  cycle  could  be         negative,  zero,  or  positive,  depending  on 
the  magnitude  of  the  applied  positive  shear  stress,  such  as  q1  or  q_  shown  on 
Fig.  12.   The  permanent  axial  strain  at  N  =  1  for  stress  reversal   '  tests  was 
always  less  than  that  developed  at  N  =  1  for  the  staged  tests,  which  had  no 
stress  reversal. 

The  rate  of  strain  development  with  number  of  cycles  was  very  clearly 
affected  by  the  existence  of  shear  stress  reversals.   Figure  13  shows  compari- 
sons of  tests  having  similar  effective  confining  pressures  and  maximum  posi- 
tive shear  stress,  both  with  and  without  shear  stress  reversal.   As  previously 
explained,  the  strains  at  N  =  1  were  always  smaller  for  the  stress  reversal 
tests.   However,  the  rates  of  strain  accumulation  were  very  much  larger  for 
the  stress  reversal  series,  although  both  types  of  test  had  positive  rates  of 
strain  accumulation. 

It  was  shown  previously  that  the  slopes  of  the  strain  versus  log  N  curves 
for  repeated  loading  without  shear  stress  reversal  were  proportional  to  the 
strain  at  the  first  cycle,  £  ,  through  the  constant,  C,  and  that  C  was  unchan- 
ged for  all  stress  levels.    This  was  not  the  case  for  the  stress  reversal 
tests.   The  slopes  of  the  e  versus  log  N  curves  for  the  stress  reversal  tests 
were  affected  by  both  the    ratio  of  cyclic  shear  stress,  Aq,  to  maximum  com- 
pressive shear  stress,  q    ,  and  the  stress  ratio  (q   I  a.) ,  where  q,.  is  the 
r    .  .       nmax       ,  .    ,  max  f  f 

static  compressive  shear   strength  for  the 

prescribed  effective  confining  pressure.   The  variation  of  the  slopes  of  the 
curves  shown  in  Fig.  13  as  a  function  of  these  stress  parameters  are  shown  in 
Fig.  14.   Note  that  Aq/q    equal  to  1.0  represents  the  condition  of  full  un- 
loading without  stress      reversal  discussed  in  the  previous  section.   Fig- 
ure 14  shows  that  as  the  compressive  failure  stress  ratios  and  amounts  of 
stress  reversal  increase,  the  slopes  of  the  e  versus  log  N  curves  increase 
rapidly.   The  trends  shown  in  Fig.  14  also     support  the  observation  from 
the  staged  tests  that  the  strain  accumulation  would  decrease  rapidly  for  par- 
tial unloadings  where  Aq/q    is  less  than  unity.   This  trend,  although  not 
completely  quantified,  was   also  observed  by  Shenton  [27]. 
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Figure  12.   Representation  of  Stress-Strain  Curves  for  Shear  Stress  Reversal 
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Figure  13.   Comparison  of  Strain  Development  Rates  for  Tests  With  and 
Without  Stress  Reversal 
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3.6   Ballast  Box  Tests 


General  Description 


In  the  laboratory  box  tests,  the  ballast  was  loaded  vertically  using  a 
simulated  section  of  tie.   The  test  setup  is  illustrated  in  Fig.  15.   The 
main  objectives  of  these  tests  were  to  investigate  the  settlement  characteris- 
tics of  the  ballast  under  field  loading  conditions,  and  to  evaluate  the  hori- 
zontal stresses  that  develop  beneath  the  tie  segment.   Details  regarding  the 
box  construction,  base  support  conditions,  instrumentation  and  complete  test 
results  can  be  found  in  Ref.  [31].   Only  the  results  which  are  pertinent  to 
analysis  in  this  report  will  be  given  here. 

The  ballast  box  was  constructed  to  represent  a  section  of  tie  about  11  in, 
long  located  directly  below  a  rail.   The  box  was  24  in.  long  by  12  in.  wide. 
The  length  was  intended  to  represent  one-half  of  the  crib  width  on  either  side 
of  the  tie  segment.   The  ballast  depth  was  12  in.  below  the  base  of  the  tie 
segment.   A  flexible  bottom  to  simulate  subgrade  support  was  incorporated  so 
that  the  deflection  basin  below  the  ballast  layer  could  be  represented.   The 
box  side  and  end  panels  were  instrumented  so  that  the  lateral  stresses  in  the 
ballast  could  be  estimated  under  both  loaded  and  unloaded  conditions.   The 
panel  numbers  and  locations  are  shown  on  Fig.  15.   The  physical  state  of  the 
ballast  prior  to  the  cyclic  loadings  was  representative  of  the  conditions 
found  under  the  tie  rail  seat  after  a  maintenance  tamping  operation. 

The  load  applied  to  the  11.5-in.  by  9-in.  tie  segment  was  4000  lb.   The 
stress  at  the  base  of  the  tie  segment  from  this  load  was  equivalent  to  the 
stress  at  the  base  of  the  tie  beneath  the  rail  seat  due  to  a  33-kip  wheel 
load. 

Permanent  Settlement 

The  permanent  settlement  trends  for  these  tests  were  similar  to  those  ob- 
served in  the  repeated  load  triaxial  tests  in  that  the  settlement  was  linearly 
related  to  the  log  of  the  number  of  cycles.   The  form  of  the  expression  for 
the  settlement  of  the  tie  for  these  ballast  box  experiments  was  identical  to 
that  developed  for  the  repeated  load  triaxial  tests  (Eq.  16),  except  that  d 
and  d  were  substituted  for  £..  and  z    ,   where  d..  and  d  were  the  residual 
vertical  deflections  of  the  tie  after  the  1st  and 

Ntn  cycles,  respectively.   Tests  on  both  compacted  and  uncompacted  samples  in- 
dicated that  the  rates  of  settlement  for  these  two  physical  states  were  differ- 
ent, but  still  both  could  be  expressed  by  the  relationship 

d^  =  dx(l  +  Cb  log  N)    ,  (17) 

where  C,  is  the  settlement  coefficient  for  the  box  tests  derived  the  same  way 
as  for   the  triaxial  tests.   The  coefficients  determined  for  the  compacted  and 
uncompacted  ballast  samples  were 

Compacted  Uncompacted 

CL   =  0.35  C,   =  0.63 

be  bu 

The  coefficient  C,   for  the  compacted  samples  was  in  general  agreement  with 
the  results  found   for  the  compacted  ballast  triaxial  tests  done  on  the  reve- 
nue site  and  FAST  granite  ballasts.   The  uncompacted  coefficient,  C   ,  was 
about  double  the  compacted  coefficient. 
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Figure  15.   Box  Test  Apparatus 
57    . 


The  physical  state  of  the  ballast  was  assessed  at  the  end  of  these  cyclic 
box  tests.   Both  ballast  density  and  plate  load  tests  showed  increases  compar- 
able to  those  measured  from  the  field  tests.   The  ballast  densities  increased 
by  about  8%  in  the  box  tests.   The  changes  between  the  pre-  and  post-mainte- 
nance ballast  densities  for  the  under-tie  locations  at  the  field  sites  were 
about  9  to  12%.   Plate  bearing  resistances  at  0.05-in.  deformation  were  com- 
pared for  the  box  tests  and  the  field  sites.   The  field  sites  showed  changes 
in  the  range  of  400  to  1200%  due  to  traffic  effects.   The  increases  in  plate 
load  resistance  measured  in  the  box  tests  were  on  the  order  of  800%  due  to  the 
cyclic  loading.   The  changes  in  ballast  density  and  bearing  stiffness  measured 
in  the  box  tests  indicated  that  the  field  loading  conditions  were  being  approx- 
imately simulated,  and  the  field  repeated-load  mechanisms  of  ballast  compac- 
tion were  being  reproduced  in  these  laboratory  tests. 

Staged  loading  tests  were  also  done  using  the  ballast  box  device.   Loads 
of  2000,  4000  and  6000  lb  were  applied  in  various  sequences.   The  4000-lb 
load  resulted  in  an  applied  stress  at  the  base  of  the  tie  segment  that  was 
equivalent  to  the  stress  at  the  base  of  a  tie  beneath  the  rail  seat  due  to  a 
33-kip  wheel  load,  as  determined  using  the  GEOTRACK  model.   The  2000-  and 
6000-lb  loads  were  selected  as  experimental  values  to  demonstrate  staged  load- 
ing effects.   The  deformation  trends  for  these  tests  were  similar  to  those  ob- 
served in  the  staged  triaxial  tests.   The  total  amount  of  vertical  deformation 
was  found  to  be  independent  of  the  applied  load  sequence.   Also,  as  in  the  tri- 
axial tests,  negligible  additional  permanent  deformation  resulted  due  to  cyc- 
ling of  loads  less  than  any  maximum  past  load. 

Horizontal  Stresses 

Estimates  of  stresses,  strains  and  deflections  for  layered  systems  such 
as  track  and  pavement  are  based  upon  an  idealized  continuum  approach  using 
either  multilayer  models,  such  as  GEOTRACK,  or  finite  element  methods.   In  a 
layered  system  where  the  stiffness  of  the  upper  layer  is  greater  than  that 
of  an  underlying  layer,  all  methods  of  analysis  predict  that  the  lower  por- 
tions of  the  stiffer  layer  will  develop  significant  incremental  tensile  stres- 
ses and  strains.   When  the  upper  layer  consists  of  an  unbound  granular  layer, 
such  as  ballast,  either  a  redistribution  of  these  stresses  must  occur  or  the 
initial  and  therefore  final  stress  state  must  be  different  than  that  predic- 
ted.  These  conclusions  are  physically  rational  since  the  unbound  layer  can- 
not sustain  tensile  stresses.   The  occurrence  of  tensile  stresses  in  an  un- 
bound granular  layer  would  result  in  very  large  inelastic  deformations  with 
each  load  application.   Such  large  inelastic  deformations  under  each  load  ap- 
plication are  not  observed  in  a  stable,  layered  system  such  as  track. 

One  of  the  primary  purposes  of  the  box  tests  was  to  measure  the  lateral 
stresses  developed  in  a  simplified  layered  system,  with  particular  interest 
in  the  development  of  residual  lateral  stresses.   It  has  been  recognized  that 
some  residual  lateral  stress  can  develop  in  soils  such  as  sands  [32]  and  clays 
[33]  due  to  compaction  and  previous  load  history,  but  the  magnitudes  of  these 
residual  stresses  have  not  been  fully  described  for  such  materials  as  ballast. 
Thus,  the  side  and  end  panels  of  the  ballast  box  were  instrumented  to  provide 
a  direct  measure  of  the  lateral  stresses  that  developed  under  the  simulated 
track  loading  conditions. 

Typical  measurements  obtained  from  the  side  panels  of  the  box  are  shown 
in  Fig.  16.   The  results  indicated  that  there  was  a  rapid  build-up  of  horizon- 
tal compressive  stresses  upon  initial  loading,  and  that  high  residual  stresses 
developed  after  only  one  cycle  of  loading.   Furthermore,  the  horizontal 
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Figure  16.  Horizontal  Stresses  on  Side  Panels  [31] 
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stresses  in  the  loaded  state  decreased  and  the  residual  stresses  increased  up 
to  about  100  cycles.   After  100  cycles,  the  horizontal  stresses  tended  to  sta- 
bilize and  the  unloaded  value  tended  to  converge  to  the  loaded  value.   This 
means  that  the  horizontal  stresses  tended  to  become  constant  during  a  load 
cycle.   The  maximum  horizontal  stresses  acting  on  the  side  panels  occurred  at 
about  the  mid-depth  of  the  ballast.   The  minimum  horizontal  stress  on  the  side 
panels  occurred  near  the  base  of  the  ballast  layer.   The  box  tests  showed  no 
evidence  of  tensile  stress  development  in  the  ballast  under  loaded  conditions. 
The  stress  measurements  from  the  instrumented  panels  also  showed  no  tendencies 
for  reduction  in  the  horizontal  stresses  in  the  ballast  due  to  the  applied  ver- 
tical load  on  the  simulated  layered  system. 

The  measured  horizontal  stresses  for  the  end  panels  are  shown  in  Fig.  17. 
Similar  trends  for  rapid  buildup  of  the  lateral  residual  stresses  were  observed 
at  the  end  of  the  ballast  box.   The  maximum  lateral  stresses  again  occurred 
about  6  in.  above  the  base  of  the  ballast  layer. 

The  variations  of  the  horizontal  stresses  under  load  and  under  no  applied 
vertical  load  as  a  function  of  the  height  above  the  box  base  are  shown  in  Fig. 
18.   In  general,  the  loaded  condition  produced  higher  lateral  stresses  for  the 
first  cycle.   After  about  100  cycles,  there  were  only  slight  differences  be- 
tween the  loaded  and  unloaded  horizontal  stress  states.   After  10,000  cycles, 
the  residual  stresses  were  almost  identical  to  the  stresses  developed  under 
full  vertical  load. 

The  measured  residual  horizontal  stresses  from  these  experiments  were 
used  to  calculate  values  of  K  ,  the  ratio  of  lateral  to  vertical  effective 
stress  at  rest  in  the  unloaded  state.   Theoretically,  the  maximum  re- 

sidual horizontal  stress  should  be  limited  by  the  passive  failure  condition. 
The  ratio  of_major  principal  stress  at  failure,  a1f»  to  minor  principal  stress 
at  failure,  o„.  at  the  passive  failure  "  condition  is  defined  by 


(18) 


for  a  strictly  cohesionless  material  such  as  ballast. 

The  calculated  K  values,  based  on  the  box  test  side  and  end  panel  meas- 
o 
urements,  the  unit  wexghts  of  the  ballast,  and  the  static  surcharge  due  to  the 

tie  segment  weight  are  shown  in  Fig.  19.   The  very  large  values_of  K  for  the 
upper  portion  of  the  side  zone  would  require  a  friction  angle,  $ ,  of  about 
56.5°,  which  might  appear  unreasonable.   However,  considering  the  curved  fail- 
ure envelopes,  possible  particle  interlocking  due  to  the  ballast  angularity, 
and  the  very  low  vertical  stresses  in  the  upper  ballast  zone,  a  friction  angle 
of  this  magnitude  may  be  possible.   Certainly  values  of  K  =6  would  be  possi- 
ble, since  the  required  friction  angle  would  be  only       about  45.5°.   The 
results  of  static  triaxial  tests  have  shown  that  this  value  would  not  be  out- 
side the  possible  range  of  values  for  low  effective  confining  pressures. 
These  experimental  observations  of  the  high  residual  lateral  stresses  and  cor- 
respondingly high  K  values  justified  the  use  of  such  values  in  the  subsequent 
GEOTRACK  analyses  for  predicting  track  response. 
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Figure  17.      Horizontal   Stresses  on  End  Panel    [31] 
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Figure  18.  Measured  Horizontal  Stress  Distribution  from  Ballast  Box 
Experiments  [31] 

62 


CVJ 


00 


vO 


CD 

CO 

CJ 

CJ 


CM 


c 

CD 

B 

•H 

CD 

CL 

X 

o 
pa 

e 
o 


a; 
Q 

U-l 

o 

a 
o 

•H 

4-1 

a 

§ 


*4 

4-1 

o 

c 
o 

•H 

4-> 

CO 
•H 

CO 

> 


Oi 


CD 
•H 


LO 


m 


(*uf)    asBg  3AoqB  aouerjsjQ 


63 


3.7   Subgrade  Material 

Sampling  Program 

Along  with  the  standard  penetration  (SPT)  and  cone  penetration  (CPT) 
tests  previously  described  in  Volume  1  [1],  undisturbed  soil  samples  (UDT) 
were  taken  at  locations  near  the  SPT  locations,  generally  only  a  few  ties 
away  along  the  centerline  of  the  track. 

In  addition  to  the  samples  recovered  from  the  field  sites,  subgrade  ma- 
terial was  recovered  from  the  FAST  track.   Undisturbed  samples  were  not  taken. 

Repeated-Load  Testing  for  Revenue  Sites 

Unconsolidated-undrained  repeated-load  triaxial  tests,  as  well  as  static 
triaxial  tests,  were  done  on  these  recovered  subgrade  materials.   The  details 
regarding  the  triaxial  testing  and  test  results  from  the  revenue  service  field 
sites,  as  well  as  index  and  classification  test  results,  can  be  found  in  Ref. 
[34]. 

The  stress  states  for  the  triaxial  repeated-load  tests  on  the  samples  re- 
covered from  the  revenue  service  sites  were  intended  to  match  the  total  stress 
states  estimated  at  the  sample  location  depths  under  representative  train  load- 
ings.  Three  basic  steps  were  necessary  to  define  these  stress  states.   The  ini- 
tial total  vertical  geostatic  stresses  at  required  depths  were  calculated  based 
on  the  average  soil  densities  as  determined  from  the  tube  samples.   The  coef- 
ficient of  lateral  stress,  K,  defined  as  the  ratio  of  total  horizontal  stress 
to  total  vertical  stress,  was  assumed  to  be  0.75  for  the  revenue  site  subgrade 
materials.   The  incremental  stresses  due  to  multiple-axle  train  loadings  were 
estimated  using  the  GEOTRACK  computer  model.   The  three-dimensional  stress 
states  at  various  depths  of  interest  were  converted  to  axisymmetric  constant 
confining  pressure  triaxial  stress  states  using  the  method  of  stress  invar- 
iants [22].   The  individual  stress  paths  used  for  these  triaxial  tests,  shown 
in  terms  of  p  =  (o  +a   )/2   and  q  =  (c,-cO/2,  can  be  found  in  Refs.  [34  and  35]. 

The  resilient  moduli  values  determined  from  these  tests  were  correlated 
with  the  cone  penetration  and  standard  penetration  test  data.   These  correla- 
tions were  necessary  since  it  was  not  possible  to  perform  a  series  of  triaxial 
tests  under  a  range  of  stress  conditions  for  each  of  the  identified  soil  stra- 
ta and  develop  stress-dependent  modular  relationships.   This  is  particularly 
true  since  remolding  of  the  undisturbed  samples  would  have  significantly  chan- 
ged their  properties.   The  cone  soundings  represented  a  continuous  measure  of 
the  soil  properties  that  could  be  used  to  extrapolate  the  properties  measured 
in  the  lab  to  other  depths.   Based  upon  the  correlations  made  and  the  test 
data  presented  in  Refs.  [34  and  35],  ranges  of  resilient  moduli  values  for 
the  subgrade  materials  were  determined.   These  ranges  of  resilient  moduli  are 
shown  in  Fig.  20  for  the  Leeds  wood  and  concrete  test  sections  and  in  Fig.  21 
for  the  Lorraine  and  Aberdeen  concrete  test  sections. 

Repeated-Load  Testing  for  FAST  Subgrade 

The  FAST  subgrade  was  also  tested  for  resilient  modulus  parameters  as  well 
as  for  index  properties  and  strength  parameters.   The  subgrade  at  FAST  is  gene- 
rally uniform  throughout  the  test  sections.   The  material  was  classified  as  a 
moderately  uniform  medium-to-fine  silty  sand.   The  Unified  Soil  Classification 
System  (USCS)  designation  of  the  FAST  subgrade  was  SM.   The  classification 
based  on  the  American  Association  of  State  Highway  and  Transportation  Officials 
(AASHTO)  was  A-2-4(0).   Since  undisturbed  samples  were  not  available,  remolded 
specimens  of  the  silty  sand  were  prepared  at  moisture  contents  of  9.4  to  10.0 
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percent. 

Two  methods  were  used  to  prepare  FAST  subgrade  soil  samples.   The  first 
sample  type  was  made  using  a  Harvard  miniature  compaction  mold.   These  sam- 
ples had  dry  unit  weights,  y,,  of  approximately  124  pcf  (101%  AASHTO  T-99) . 
The  other  sample  type  was     made  using  a  specially  fabricated  sample  mold 
and  using  a  preparation  technique  similar  to  that  described  by  Ladd  [36]  with 
which  the  soil  density  state  could  be  controlled.   The  resulting  dry  unit 
weights,  y , ,  using  the  controlled  compaction  technique  were  112  pcf  (91% 
AASHTO      T-99).   All  soil  samples  had  length-to-diameter  ratios  of  2.5:1. 
Specific  details  regarding  testing  procedures,  index  test  results,  and  strength 
parameters  for  the  FAST  subgrade  are  given  in  Ref.  [37]. 

The  results  from  the  FAST  subgrade  tests  indicate  that  the  type  of  compac- 
tion used  to  prepare  samples  can  have  an  effect  on  the  measured  resilient  modu- 
li values.   Samples  prepared  using  the  undercompaction  method  resulted  in  gene- 
rally higher  moduli  than  those  prepared  using  the  Harvard  miniature  method. 
However,  the  data  from  both  of  these  methods  were  combined  and  used  to  deter- 
mine the  modulus  coefficient,  K-.  ,  and  modulus  exponent,  K~,  for  the  following 
resilient  formulation: 

0  K2 
a 

A  vibratory  method  of  compaction  was  used  to  prepare  the  subgrade  samples  test- 
ed in  Ref.  [17].   Values  for  K  and  K  were  also  obtained  for  the  data  from 
this  source.   The  resulting  parameters  are  as  follows: 

Source  1  2 

Stewart  877  1.10 

Thompson  [17]  654  1.08 

The  main  difference  is  in  the  coefficient  of  K. . 

A  possible  explanation  for  this  difference  is  that  the  resilient  moduli 
values  taken  after  40,000  cycles  were  used  from  the  tests  in  this  study,  where- 
as the  values  reported  in  Ref.  [17]  were  based  on  staged  test  results  taken 
after  5000-load  repetitions  following  a  stress  state  change.   It  was  previous- 
ly shown  that  single  stage  testing  generally  results  in  increasing  modulus  with 
increasing  number  of  cycles,  but  staged  testing  does  not  show  a  continuous  in- 
crease with  number  of  cycles. 
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CHAPTER  4.   TRACK  MODULUS 


4.1   Introduction 


The  use  of  track  modulus  for  assessing  track  performance  is  common  in  the 
railroad  industry.   The  theoretical  formulation  of  track  modulus  is  based  on 
the  assumption  that  the  rail  acts  as   a  beam  continuously  supported  on  an 
elastic  foundation.   The  track  modulus  is  defined  as  the  supporting  force  per 
unit  length  of  rail  per  unit  deflection  in  the  track  system.   The  resulting 
equation  for  track  modulus  is 


*y 


pT~   T1 


'(6°   '  EI     '  (20) 

where  u  is  the  track  modulus  (units  =  F/L/L) , 

P  is  the  applied  load  (units  =  F) , 

6  is  the  vertical  rail  deflection  (units  =  L) ,  and 

EI  is  the  rail  bending  stiffness  (units  =  FL^). 

An  important  difference  between  the  actual  track  support  and  the  ideali- 
zed formulation  of  a  rail  on  an  elastic  support  is  that  the  rail  load  is  actu- 
ally applied  to  the  foundation  through  discrete  supports,  which  are  the  ties, 
rather  than  through  support  distributed  along  the  track  foundation.   Another 
difference,  for  concrete  tie  track  systems,  is  the  inclusion  of  flexible  tie 
pads  between  the  rails  and  the  tie  rail  seats.   Even  though  these  differences 
exist  between  the  actual  track  structure  and  the  theoretical  formulation,  the 
track  modulus  has  historically  been  used  as  a  measure  of  track  quality.   Fur- 
ther details  on  the  historical  development  and  interpretation  of  track  modulus 
can  be  found  in  Refs.  [38  and  39]. 

This  chapter  will  present  the  track  structure  and  foundation  properties 
selected  for  input  into  the  GEOTRACK  program  in  order  to  predict  the  track 
modulus  at  each  of  the  revenue  field  sites.   These  predictions  will  be  com- 
pared with  field  measurements,  and  the  trends  will  be  discussed. 

4.2   Selection  of  Parameters 

The  GEOTRACK  model  [40]  was  used  to  determine  elastic  track  deflections 
for  predicting  values  of  track  modulus  for  the  revenue  field  sites.   The  GEO- 
TRACK computer  program  is  a  three-dimensional,  multi-layer  model  for  determin- 
ing the  elastic  response  of  the  track  structure,  using  stress-dependent  proper- 
ties for  the  ballast,  subballast, and  subgrade  materials.   The  output  of  the 
program  provides  rail  seat  loads,  tie-ballast  reactions,  tie  and  rail  deflec- 
tions, and  tie  and  rail  bending  moments.   In  addition,  the  output  provides  ver- 
tical displacements  and  the  complete  three-dimensional  stress  state  at  selec- 
ted locations  in  the  ballast,  subballast  and  subgrade.   The  validity  of  this 
model  for  evaluating  the  resilient  response  of  track  was  established  by  compar- 
ing measured  foundation  responses  obtained  from  FAST  instrumentation  with  pre- 
dictions made  using  the  GEOTRACK  model.   Details  regarding  the  model  valida- 
tion, measured  FAST  foundation  responses,  and  the  results  of  parametric  studies 
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done  using  the  GEOTRACK  model  can  be  found  in  Refs.  [41  and  42]. 

In  order  to  make  determinations  of  the  elastic  responses  of  the  revenue 
sites,  the  structural  and  foundation  properties  were  determined  for  each  of 
the  sites.   The  structural  properties  are  given  in  Table  2.   Figures  22  to  25 
show  the  ballast  and  subgrade  layer  divisions  and  properties. 

The  shear  stress-resilient  strain  formulation  developed  in  Chapter  3  was 
used  to  characterize  the  stress-dependent  resilient  moduli  for  these  sites. 
The  other  ballast  parameter  values  are  shown  in  Figs.  22  to  25. 

The  subgrade  layer  properties  used  for  the  GEOTRACK  analyses  were  chosen 
based  on  the  results  shown  in  Figs.  20  and  21.   Layer  divisions  for  the  sub- 
grades  were  made  (Figs.  22  to  25)  where  there  appeared  to  be  significant  chan- 
ges in  the  measured  resilient  properties.   The  average  resilient  modulus  for 
each  subgrade  layer  was  used  as  the  representative  value  for  the  layer.   The 
moduli  for  the  subgrade  layers  were  held  constant,  since  stress-state-depen- 
dent relationships  were  not  available  for  the  subgrade.   The  values  of  total 
unit  weight  were  determined  from  the  subgrade  samples  recovered  from  the  un- 
disturbed tube  sampling.   The  values  of  Poisson's  ratio  used  for  the  subgrade 
were  based  upon  static  triaxial  tests  and  were  generally  very  close  to  0.5. 
The  ratio  of  lateral  to  vertical  stress,  K  =  0.75,  was  estimated  from  the  soil 
conditions  as  a  reasonable  value  for  the  subgrade  layers.   An  exception  to 
these  values  was  that  a  K  =2  value  was  used  for  the  upper  granular  subgrade 
at  the  Leeds  sites. 

The  parameter  K  is  defined  as  the  ratio  of  lateral  effective  stress  to 
vertical  effective    stress  in  the  granular  layers  under  the  condition  of 
zero  lateral  strain.   The  values  of  K  =  4  to  6  in  the  upper  granular  layers 
may  appear  to  be  unusually  high.       However,  these  values  were  based  upon 
experimental  observations  of  ballast  behavior  in  a  laboratory  ballast  box  de- 
vice developed  for  simulating  field  conditions. 

4.3   Comparison  with  Field  Measurements 

Measurements  of  track  vertical  load-deflection  response  were  made  at  each 
of  the  revenue  field  sites  by  Battelle-Columbus  Laboratories  (BCL) .   For  the 
test,  a  point  load  of  known  magnitude  was  simultaneously  applied  to  both  rails, 
using  vertical  hydraulic  jacks  reacting  against  a  loaded  freight  car.   Rail  de- 
flections were  measured  by  sighting  a  steel  scale  with  a  surveyor's  theodolite. 
The  scale  was  attached  to  the  rail.   Seven  to  ten  measurements  at  random  tie 
locations  were  made  to  assess  the  variability  in  the  load-deformation  respon- 
ses within  the  track  sections.   These  measurements  were  made  both  before  and 
after  the  maintenance  operations  to  see  what  differences  in  track  support  re- 
sulted from  the  surfacing  operation. 

Track  load-deflection  measurements  were  made  at  7  to  10  tie  locations  at 
each  site.   The  measured  averages  and  standard  deviations  for  the  pre-  and 
post-maintenance  load-deflection  curves  were  given  in  Volume  1  [1] .   A  small 
amount  of  slack  between  the  bottom  of  the  tie  and  the  top  of  the  ballast  layer 
may  have  contributed  to  the  nonlinear  response  in  the  initial  portion  of  the 
load-deflection  curves.   Most  of  this  initial  slack  was  assumed  to  have  been 
eliminated  after  about  a  6-kip  load  was  applied.   To  minimize  this  effect,  the 
track  moduli  values  used  for  this  comparison  were  calculated  for  the  6-  to  30- 
kip  load  range. 

The  mean  track  moduli  measurements  and  estimates  of  the  standard  devia- 
tion were  made  using  the  mean  load-deflection  curves  and  the  deflections  at 
±1  standard  deviation,  both  taken  at  the  6-  and  30-kip  load  levels.   The 
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standard  deviations  of  track  modulus  calculated  in  this  way  are  not  symmetri- 
cal about  the  mean  value,  and  slightly  underestimate  the  actual  standard  devia- 
tions.  For  the  purposes  addressed  in  this  chapter,  the  standard  deviations  of  the 
track  moduli  values  will  refer  to  the  limits  calculated,  based  on  the  standard 
deviations  of  the  rail  deflections. 

To  determine  track  modulus  with  GEOTRACK,  deflections  were  calculated  for 
the  single  axle  solution  with  loads  of  6  and  30  kips.   The  difference  in  loads 
and  difference  in  deflections  were  substituted  into  Eq.  20  to  get  track  modu- 
lus. 

The  measured  values  of  track  modulus  and  the  estimated  standard  devia- 
tions for  all  of  the  revenue  field  sites  are  given  in  Table  16  and  shown  in 
Fig.  26,  along  with  the  predicted  values  based  on  the  GEOTRACK  analyses.   Seve- 
ral items  on  Fig.  26  deserve  attention.   First,  there  was  no  significant  change 
in  the  track  moduli  values  as  a  result  of  the  surfacing  operations.   However, 
the  pre-maintenance  values  did  seem  to  be  less  variable  than  the  post-mainte- 
nance values.   This  variation  was  also  apparent  in  the  average  load-deflection 
curves,  where  the  scatter  about  the  mean  was  visibly  larger  for  the  post-main- 
tenance values.   The  fact  that  the  average  measured  values  were  greater  in  the 
Leeds  wood  and  Lorraine  sections  after  maintenance  is  probably  not  statisti- 
cally significant,  since  the  estimated  standard  deviations  all  overlap. 

Another  observation  from  Fig.  26  is  that  there  did  not  appear  to  be  a  cor- 
relation between  height  of  raise  and  post-maintenance  track  modulus.   Raises  of 
1.5  to  2  in.  were  given  to  the  Leeds  and  Lorraine  sections  and  only  about  0.1 
in.  to  the  Aberdeen  section.   In  spite  of  this,  the  pre-  and  post-maintenance 
values  for  any  one  test  section  were  approximately  equal,  and  the  Aberdeen  val- 
ue was  between  the  Lorraine  and  Leeds  values. 

An  explanation  for  the  increased  variability  of  the  post-maintenance  modu- 
li values  as  compared  to  the  pre-maintenance  values  could  be  that  the  surfacing 
decreased  the  uniformity  of  track  support  conditions  between  tie  locations. 
One  purpose  of  track  maintenance  is  to  improve  the  overall  track  surface,  mean- 
ing to  smooth  out  the  vertical  track  roughness.   A  larger  amount  of  raise  must 
be  applied  beneath  the  lowest  ties  to  achieve  a  uniform  surface.   These  varia- 
tions in  the  actual  raises  applied  beneath  the  individual  ties  could  cause  lo- 
cal differences  in  the  amount  of  ballast  disturbance,  hence  variations  in  the 
ballast  physical  state  from  one  tie  to  another. 

Part  of  the  differences  in  absolute  magnitude  of  the  average  field  track 
modulus  measurements  can  be  explained  in  terms  of  the  differences  in  the  track 
substructures.   The  parametric  study  [41]  using  .the  GEOTRACK  model  indicated 
that  track  modulus  increased  as  ballast  depth  increased.   The  ballast  profiles 
and  the  simplified  layer  characterizations  used  in  the  GEOTRACK  analyses  (Figs. 
22  to  25)  show  that  the  Leeds  wood  section  had  only  about  9  in.  of  ballast  be- 
neath the  tie,  whereas  the  Leeds  concrete  section  contained  about  14  in.   The 
Aberdeen  site  had  20  in.  of  ballast,  and  the  ballast  depth  at  Lorraine  was  es- 
timated to  be  28  in.  below  the  tie.   This  trend  of  increasing  track  modulus 
with  increasing  ballast  depth  for  the  field  sites  is  clearly  confirmed  by  the 
field  measurements  in  Fig.  27. 

The  track  modulus  is  also  influenced  by  the  subgrade  characteristics. 
The  GEOTRACK  model  indicates  that  the  compression  of  the  ballast  layer  accounts 
for  about  10  to  20%  of  the  total  vertical  rail  deflection.  Most  of  the  remain- 
ing 80  to  90%  of  the  rail  deflection  is  accounted  for  by  the  compression  of  the 
subgrade  materials.  Furthermore,  25  to  40%  of  the  subgrade  deformation  indica- 
ted by  GEOTRACK  occurs  below  a  depth  of  about  10  ft,  even  though  the  stresses 
below  this  depth  are  low  [41] . 
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The  Lorraine  test  section  was  found  to  have  the  greatest  depth  of  ballast- 
type  material  and  the  stiff est  subgrade.   Correspondingly,  the  Lorraine  sec- 
tion had  the  highest  values  of  measured  and  predicted  track  modulus  of  these 
four  revenue  sites.   As  can  be  seen  in  Fig.  27,  the  predicted  value  of  track 
modulus  was  higher  than  the  average  measured  values,  but  well  within  the  esti- 
mated standard  deviations. 

The  Aberdeen  post-maintenance  value  was  lower  than  the  Lorraine  value  due 
to  a  combination  of  reduced  ballast  thickness  and  lower  overall  subgrade  stiff- 
ness.  The  predicted  track  modulus  for  the  Aberdeen  site  was  in  very  close 
agreement  with  the  measured  values. 

For  the  Leeds  sites,  the  predicted  values  of  track  modulus  were  higher 
than  the  field  values.   However,  the  field  values  appear  to  be  unusually  low. 
The  lower  ballast  thickness  at  the  Leeds  sites  can  account  for  some  of  the 
difference  between  the  Leeds  sites  and  the  other  two  field  sites,  but  these 
differences  in  ballast  layer  thickness  were  not  enough  to  cause  the  low  values 
measured  at  both  Leeds  sites. 

Considering  the  similarity  between  the  subgrade  stiffnesses  at  the  Leeds 
sites  and  the  Aberdeen  site,  closer  agreement  between  the  field  measurements 
from  these  sites  would  be  expected.   It  is  possible  that  the  embankment  in  the 
Leeds  wood  section  resulted  in  reduced  subgrade  confinement  and  hence  increas- 
ed vertical  deflections.   This  would  result  in  lower  values  of  track  modulus. 
However,  the  concrete  section  at  Leeds,  which  was  built  at  grade,  had  compara- 
bly low  track  moduli  values  as  well.   Thus,  the  embankment  condition  must  not 
have  been  a  major  factor. 

It  would  be  necessary  to  reduce  the  subgrade  stiffnesses  that  were  used 
in  the  GEOTRACK  analyses  by  at  least  50%  in  order  to  match  the  field  track 
modulus  measurements  at  the  Leeds  test  sections.   However,  no  rational  justi- 
fication for  making  adjustments  of  these  magnitudes  could  be  found.   If  the 
soil  moduli  were  overestimated  at  the  Leeds  sites,  a  similar  systematic  error 
should  have  occurred  with  the  Lorraine  and  Aberdeen  subgrades.   Since  the  pre- 
dicted Lorraine  and  Aberdeen  values  were  in  good  agreement  with  the  measure- 
ments, a  similar  adjustment  to  the  subgrade  moduli  at  those  sites  would  shift 
the  predicted  values  away  from  the  measurements. 

The  measured  and  predicted  values  of  track  modulus  are  dependent  on  seve- 
ral factors,  one  of  which  is  the  support  condition  of  the  tie.   These  support 
conditions  are  a  function  of  track  settlement  and  maintenance  effects.   The 
support  conditions  would  affect  concrete  and  wood  ties  in  a  different  fashion. 

The  GEOTRACK  model  uses  uniform  properties,  throughout  each  layer,  includ- 
ing the  support  and  contact  under  all  ties.   The  field  plate  load  tests  showed 
that  the  ballast  stiffness  was  not  uniform  under  the  ties,  but  was  greater 
near  the  rail  seat  areas  than  under  the  tie  centers.   Furthermore,  after  the 
maintenance  the  ballast  physical  state  was  more  uniform  than  before  mainte- 
nance.  However,  even  uniform  physical  state  or  modulus  does  not  result  in 
uniform  support  conditions  along  the  tie. 

After  a  maintenance  operation  in  which  a  high  raise  was  given  to  the 
track,  the  ballast  physical  properties  would  be  more  uniform.   However,  the 
tamping  and  raise  during  the  maintenance  is  done  only  near  the  rail  seats, 
which  could  cause  a  gap  near  the  center  of  the  ties.   This  causes  the  actual 
load  bearing  areas  to  be  near  the  rail  seats  due  to  the  lack  of  contact  near 
the  tie  centers.   This  increased  load  bearing  near  the  rail  seats  would  cause 
higher  rail  deflections  directly  under  the  applied  load  than  would  result  from 
a  continuously  supported  condition.   The  GEOTRACK  model  uses  continuous  con- 
tact between  the  tie  directly  under  the  applied  vertical  load  and  the  ballast 
surface,  leading  to  a  lower  rail  deflection  and  higher  track  modulus  than 
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measured  in  the  field  after  maintenance. 

As  traffic  accumulates  over  the  track,  the  ballast  beneath  the  rail  seats 
is  recompacted  and  becomes  stiffer,  but  the  tie  contact  approaches  a  continuous 
condition.   The  more  uniform  support  across  the  tie  leads  to  a  lower  rail  de- 
flection and  hence  a  higher  track  modulus.   After  a  small  maintenance  raise, 
the  initial  ballast  physical  states  may  be  less  uniform  under  the  tie  than 
after  a  large  raise.   However,  in  the  former  case,  the  contact  may  be  more 
continuous  and  the  ballast  will  have  greater  stiffness.   This  would  result  in 
a  higher  track  modulus  after  a  small  raise  than  after  a  large  raise,  or  at 
least  not  much  change  in  the  pre-  and  post-maintenance  values. 

Further  traffic  leads  to  a  centerbound  track,  in  which  a  greater  portion 
of  the  rail  load  would  be  carried  by  the  central  portion  of  the  ties  than  in 
non-centerbound  track.   A  gap  also  tends  to  develop  between  the  ballast  and 
tie  beneath  the  rail.   If  the  gap  is  included  in  the  deflection  measurement, 
then  the  calculated  track  modulus  would  appear  to  decrease  as  a  result  of  cen- 
terbinding.   However,  if  the  gap  is  eliminated  by  a  seating  load,  then  the 
track  modulus  would  probably  increase  with  centerbinding. 

The  interactions  between  the  variable  ballast  physical  states,  tie  sup- 
port conditions,  and  structural  factors  such  as  tie  stiffness  and  rail  size 
make  generalizations  about  track  modulus  uncertain.   This  is  particularly  true 
since  the  degree  and  type  of  maintenance  disturbance  and  traffic  history  of  a 
site  can  change  the  ballast  physical  states  in  varying  amounts.   The  scatter 
of  the  field  measurements  were  such  that  there  were  no  clear  trends  distin- 
guishing the  pre-maintenance  track  moduli  values  from  the  post-maintenance 
values.   The  predictions  of  track  modulus  using  the  GEOTRACK  program  are  some- 
what limited  by  the  uniform  layer  property  and  full  contact  representations, 
and  the  inability  to  represent  the  maintenance  factors  for  the  field  sites. 
For  these  reasons,  variations  between  the  measured  and  predicted  track  moduli 
values  for  the  sites  can  be  expected  due  to  the  variations  in  ballast  proper- 
ties and  support  conditions  which  were  affected  by  the  maintenance  operations. 

The  possible  centerbinding  and  uniformity  of  support  conditions  beneath 
the  tie  bottom  may  not,  however,  be  a  significant  contributing  factor  to  the 
track  modulus  values,  although  the  effects  are  physically  rational.   Differen- 
ces between  the  bending  stiffnesses  of  the  wood  and  concrete  ties  would  also 
not  result  in  large  differences  in  track  moduli.   Since  only  10  to  20%  of  the 
total  track  deformation  occurs  due  to  ballast  compression,  the  subgrade  deflec- 
tions appear  to  be  much  more  important.   It  was  shown  in  previous  parametric 
studies  [41,42]  that  variations  in  tie  stiffnesses  do  not  have  a  major  effect 
on  the  vertical  subgrade  stresses  beneath  the  rail  seats.   The  effects  of  a 
concrete  tie  (EI  -  1,740,000  to  2,360,000  K-in.2)  have  shown  that  the  vertical 
stress  beneath  the  rail  seat  may  be  only  about  1  psi  greater  than  the  wood  tie 
(EI  -  386,000  K-in. 2)  for  similar  conditions.   Thus,  the  subgrade  contribution 
to  track  modulus  should  be  about  the  same  for  wood  and  concrete  ties.   The  FAST 
dynamic  measurements  of  resilient  subgrade  deflection  also  showed  no  signifi- 
cant difference  between  the  deflections  in  wood  and  concrete  sections.   This 
would  indicate  that  although  the  stress  distributions  and  deformations  in  the 
ballast  layer  were  affected  by  tie  stiffness  and  possibly  centerbinding,  the 
subgrade  responses  were  controlled  mainly  by  the  subgrade  properties  for  the 
same  ballast  layer  thickness. 

A  comparison  between  the  measured  and  predicted  values  of  the  revenue 
site  vertical  rail  deflection  is  shown  in  Fig.  28.   The  differential  deflection 
equals  the  deflection  at  30  kips  wheel  load  minus  the  deflection  at  6  kips.   A 
change  in  rail  deflection  causes  a  larger  change  in  the  calculated  value  of 

79 


0.12  r- 


0.10  - 


•H 

o  0.08 

r-t 

Q 


•h  0.06  " 


T3 

0) 
H 
3 
CO 
(0 

2 


0.04  - 


0.02 


0.02       0.04       0.06       0.08 
Predicted  Differential  Rail  Deflection  (in.) 


0.10 
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track  modulus.   For  example,  the  predicted  values  of  track  modulus  at  the 
Leeds  wood  and  concrete  sections  were  in  error  by  about  60  to  100%,  but  the 
errors  in  predicted  rail  deflections  were  less.   Furthermore,  the  differences 
between  predicted  and  measured  deflections  were  small  in  absolute  magnitude. 
Thus,  the  predicted  values  of  the  rail  deflections  were  in  reasonable  agree- 
ment with  the  measured  values.   In  fact,  considering  the  complexity  of  soil 
behavior,  the  accuracy  of  deformation  prediction  obtained  was  very  good. 
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CHAPTER  5.   TRACK  LOADING 


5.1   Introduction 


Any  methodology  for  predicting  the  accumulated  permanent  strains  and  set- 
tlement of  track  must  take  into  account  the  actual  loading  conditions.   The 
actual  loading  of  a  given  track  will  consist  of  a  mix  of  many  different  wheel 
loads  which  are  defined  by  individual  car  weights  and  the  dynamic  response  of 
the  moving  train. 

The  approach  taken  in  this  study  for  determination  of  the  mix  of  wheel 
loadings  at  the  sites  was  based  on  the  direct  measurement  of  the  actual  wheel 
loads  at  each  site,  rather  than  the  use  of  empirical  load  factors.   The  fol- 
lowing sections  will  present  these  measurements,  and  describe  the  method  used 
to  characterize  the  mix  of  wheel  loads  present  at  each  of  the  field  sites,  in- 
cluding FAST.   The  static  wheel  load  distributions  for  each  of  the  sites  were 
obtained  from  the  railroad  companies,  based  on  the  car  weights  and  tonnage 
records. 

5.2  Measured  Wheel  Loads 

Strain  gauge  instrumentation  was  installed  on  the  rails  to  measure  the 
actual  wheel  loads.   This  work  was  performed  by  BCL.   Descriptions  of  the 
strain  gauges,  the  calibration  procedures,  the  instrument  locations,  and  the 
data  recording  methods  are  given  in  [44]  and  in  Volume  3  [43].   On  the  basis  of 

past  experience  in  wheel  load  detection,  a  sufficient  number  of  loads  were 
recorded  at  each  site  so  that  the  load  distributions  were  accurately  defined. 
The  data  were  reduced  to  the  form  of  the  percent  of  wheel  loads  that  ex- 
ceeded a  given  load  level.   This  is  analogous  to  the  cumulative  distribution 
function  (CDF)  of  random  variables.   The  random  variables  in  this  case  were 
the  static  and  dynamic  wheel  loads.   The  percent  exceeded  level  for  the  meas- 
ured wheel  loads  was  converted  to  the  CDF  by 

Fx(x)  =  1QQ  -  P^t  exceeded  =  p(x  <  x)    ,  (n) 

where  Fy(x)  is  the  probability  that  the  given  wheel  load,  X,  is  assigned  a 
value    equal  to  or  less  than  a  specified  value,  x.   This  statement  is  equi- 
valent to  the  expression  P(X  <^  x)  . 

The  cumulative  distributions  for  the  static  and  dynamic  wheel  loads  were 
given  in  Volume  1  [1]  for  the  field  sites  at  Leeds,  Lorraine,  Aberdeen, and 
FAST.   Several  trends  were  identified  in  these  distributions.   There  were  sig- 
nificant differences  between  the  low  frequency  of  occurrence  wheel  loads  at 
the  sites.   These  differences  appeared  as  wide  separations  between  the  static 
and  dynamic  CDF's  near  the  upper  tails  of  the  distributions.   These  separa- 
tions began  at  about  the  80-90  percent  exceeded  level  (CDF  =  0.8  -  0.9)  for 
the  revenue  sites.   The  largest  variation  between  the  low  frequency  of  occur- 
rence static  and  dynamic  wheel  loads  was  at  the  Aberdeen  site.   The  very  high 
dynamic  loads,  a  factor  of  3  to  5  times  greater  than  the  mean  dynamic  load, 
are  thought  to  be  caused  by  tread  surface  irregularities  and  the  high  speeds 
of  the  passenger  trains. 
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The  wheel  load  distributions  for  the  Lorraine  site  showed  two  distinct 
traffic  patterns.   The  breaks  in  the  wheel  load  distribution  curves  for  this 
site  were  due  to  the  fact  that  approximately  equal  numbers  of  loaded  and  empty 
trains  pass  over  this  site.   Most  of  this  traffic  consisted  of  unit  coal 
trains  that  traveled  full  in  one  direction  and  returned  empty  along  the  single 
mainline  track. 

The  differences  between  the  measured  static  and  dynamic  wheel  loads  were 
not  as  severe  at  FAST.   The  relatively  close  agreement  between  the  static  and 
dynamic  wheel  loads  at  FAST  was  due  to  the  test  control  and  well  maintained 
equipment  at  that  site.   For  example,  tread  surface  irregularities  are  regular- 
ly repaired  at  FAST,  eliminating  a  major  cause  of  high  dynamic  impact  loads. 
In  addition,  the  static  wheel  loads  at  FAST  were  very  uniform  for  the  consist 
used. 

On  the  basis  of  the  data  collected  by  BCL,  the  mean  dynamic  and  static 
wheel  loads  for  each  of  the  sites  were  calculated  according  to  methods  des- 
cribed in  Ref.  [45].   These  mean  static  and  dynamic  wheel  loads  are  given  in 
Table  17  for  the  revenue  sites  and  FAST. 

5.3  Load  Distribution  Methodology 

Selection  of  Wheel  Loads 

Although  the  complete  dynamic  load  distributions  were  defined  for  the 
sites,  it  would  not  be  practical  to  determine  the  effect  that  the  continuous 
load  spectrum  would  have  on  the  track  structures.   To  simplify  the  continuous 
load  distributions  and  reduce  the  number  of  computer  analyses  necessary,  dyna- 
mic wheel  loads  corresponding  to  the  50,  10,  1  and  0.1  percent  exceeded  levels 
were  chosen  to  represent  the  entire  dynamic  load  distributions.   The  meaning 
of  the  0.1  percent  exceeded  level,  for  example,  is  that  only  1  out  of  1000 
loads  will  be  greater  than  the  specified  value.   The  CDF  corresponding  to  the 
50,  10,  1  and  0.1  percent  exceeded  levels  were  0.5,  0.9,  0.99  and  0.999,  res- 
pectively.  More  rigorous  methods  have  been  suggested  to  represent  the  loadings 
for  highway  systems,  but  the  chosen  levels  of  exceedance  are  sufficient  to  ade- 
quately cover  the  range  of  expected  dynamic  wheel  loads. 

Once  these  CDF  values  were  chosen,  it  was  necessary  to  determine  the  rela- 
tive frequencies  of  occurrence  for  these  dynamic  load  levels.   The  method  used 
was  such  that  the  midpoint  loads  between  two  values  of  the  CDF  were  those  loads 
corresponding  to  the  50,  10,  1  and  0.1  load  levels.   The  CDF  ranges  for  each 
load  level  were  chosen  such  that  the  average  CDF  values  corresponded  to  the  se- 
lected values  of  0.5,  0.9,  0.99  and  0.999.   Thus,  the  difference  between  any 
two  successive  values  of  the  CDF  was  used  to  represent  the  relative  frequency 
of  occurrence  for  the  given  dynamic  load.   The  procedure  used  was  as  follows: 

1)  Begin  at  the  upper  limit  of  the  CDF,  where  F  (°°)  =  1. 

2)  Calculate  the  next  value  of  F  (x- )  such  that 

A   -L 

F  (oo)  -  F  (x  ) 

— 2  =     °*999    *  (22) 

Note  that  Fx(oo)  -  Fx(x  )  =  ?(x±  <_  X±   £  «>)  =  f^*-^ 

where  f  (x  )  =  relative  frequency  of  occurrence  for  the  0.1  percent  ex- 
ceeded load,  and  X..  =  dynamic  wheel  load  at  the  0.1  percent  exceeded 
level. 
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3)  Calculate  the  next  value  of  F  (x  )  such  that 

F  (x  )  -  F  (x  ) 

-^— - 2  =   °'"    '  (23) 

Note  that  Fx(x  )  -  Fx(x2)  =  P(x2<  X2  <_   x  )  -  fx(x2) 

where  f  (x~)  and  X  are  the  relative  frequency  of  occurrence  and  dyna- 
mic wheel  load  for  the  1  percent  exceeded  level. 

4)  Continue  this  process  for  CDF  =  0.9  and  0.5. 

5)  Since 

1 

the  remainder  of  the  loads  below  the  50  percent  exceeded  level  must  be 

accounted  for.   To  do  this,  calculate  f  (x  )  such  that 

X   _) 

4 

fx(x5)  =  1-1   f(x.)    ,  (25) 

i=l 

where  f  (x  )  =  relative  frequency  of  loads  less  than  the  50  percent  ex- 
x  j 

ceeded  level.   The  dynamic  load,  X,.»  is  the  load  where  the  CDF  has  a 

value  given  by 

F  (x  )  -  F  (-oo) 
Fx(x5)  =  -^—2 -± =  0.091    ,  (26) 

where  F  (-oo)  =  0. 
x 

This  method  is  somewhat  arbitrary,  but  the  approach  is  rational  and  should 

approximately  represent  the  mix  of  wheel  loads.   The  dynamic  loads  for  the  0.1, 

1,  10  and  50  percent  exceeded  levels  for  the  field  sites  and  FAST,  along  with 

the  relative  frequencies  calculated  using  the  above  outlined  procedure,  are 

given  in  Table  18. 

Number  of  Load  Applications 

In  order  to  relate  laboratory  repeated  load  tests  to  the  traffic  levels  at 
the  sites,  it  was  necessary  to  convert  MGT  values  to  an  equivalent  number  of 
cycles  of  load  application.   The  dynamic  records  obtained  from  instrumented  sec- 
tions of  the  FAST  track  were  studied  [46],  and  it  was  concluded  that  each  axle 
pair  in  a  truck  acts  nearly  as  a  single  load  cycle  for  the  ballast  layer  of  the 
track  support  system.   At  greater  depths,  the  stress  distributions  from  the  ad- 
jacent axles  overlap  sufficiently  so  that  the  four  axles  from  adjacent  wheel 
trucks  represent  a  single  load  cycle.   Therefore,  two  axles  will  be  used  to  rep- 
resent a  load  cycle  for  the  ballast  layers  and  four  axles  to  represent  a  single 
load  cycle  for  the  subgrade  layers. 

Now  1  MGT  equals  2  x  10^  kips,  and  the  number  of  kips  per  axle  equals 
twice  the  wheel  load  in  kips.   Thus  the  number  of  axles  per  MGT  is 

106 
N  =  t^t   >  (27) 


a   E(X) 
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where  N  =  number  of  axles  per  MGT,  and 
a 

E(X)  =  expected  or  mean  value  of  the  static  wheel  load  distribution,  in 
kips. 

The  mean  static  wheel  load  (Table  17)  was  used  rather  than  the  mean  dynamic 
load,  since  railroad  tonnage  records  are  based  upon  the  actual  weights  of  the 
cars. 

The  train's  axle  load  in  tons  is  numerically  equal  to  the  wheel  load  in 
kips,  assuming  equal  distribution  of  load  on  both  wheels.   Since  the  axle  load 
acts  on  both  rails,  the  number  of  wheels  passing  a  particular  location  per  MGT 
is  the  same  as  that  calculated  using  Eq.  27. 

Based  upon  the  previous  conclusions  regarding  the  number  of  axles  per  cy- 
cle for  the  ballast  and  subgrade  layers,  the  conversion  from  MGT  to  number  of 
load  applications  for  the  ballast  is 

N 
\  =  y"       ,  (28) 

and  for  the  subgrade  is 

N 

s    4 

where  N,  and  N  are  the  number  of  load  applications  per  MGT  for  the  ballast  and 
subgrade  layers.   Using  the  mean  static  wheel  load  for  the  sites  as  given  in 
Table  17,  the  equivalent  numbers  of  load  applications  per  MGT  for  the  indivi- 
dual field  sites  were  calculated.   These  results  are  summarized  in  Table  19. 

The  distributions  of  the  dynamic  loads  as  previously  calculated  were  used 
to  define  the  mix  of  wheel  loads  at  each  of  the  sites.   To  determine  the  num- 
ber of  axle  loads  of  the  prescribed  magnitudes,  the  relative  frequency  for 
that  load  level,  f  (x.),  is  multiplied  by  the  equivalent  number  of  cycles  per 
MGT,  as  given  by 

nm  =  WNb     •  (30) 

for  the  ballast,  and  by 

Nsl=£x(x.)Ns    ,  (31) 

for  the  subgrade.   The  equivalent  number  of  load  cycles  per  MGT,  at  the  pre- 
scribed dynamic  wheel  loads,  is   given  in  Table  20  for  both  the  ballast  and 
subgrade  layers.   Again,  note  that  the  single  wheel  load  in  kips  is  equal  to 
the  axle  load  in  tons. 

Multiple-Axle  Loadings 

The  GEOTRACK  computer  model,  to  be  used  later  in  analyzing  the  track  per- 
formance, was  modified  during  this  program  to  handle  multiple-axle  loads.   A 
four-axle  loading  condition,  representing  two  adjacent  wheel  trucks,  was  used 
to  represent  the  dynamic  loadings  at  the  field  sites.   Since  it  was  unlikely 
that  actual  train  loadings  would  consist  of  several  wheel  loads  of  the  same 
magnitude  occurring  sequentially,  the  complete  representation  of  dynamic  track 
loadings  was  defined  by  combining  three  axles  of  the  mean  dynamic  load  with  one 
axle  from  the  individual  values  of  representative  dynamic  loads. 
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Table  19.   Number  of  Cycles  per  MGT  for  Ballast  and  Subgrade 

Leeds    Lorraine   Aberdeen     FAST 

Ballast,  Nu      31000      28000      27000      17000 
b 

Subgrade,  N      16000      14000      13000       8000 


Table  20.   Number  of  Cycles  per  MGT  at  Prescribed  Dynamic  Wheel  Loads 
for  Ballast  and  Subgrade 


P  (kips) 

6.2 

15.0 

27.3 

35.7 

43.5 

Leeds 

Nb 

5700 

20000 

5200 

500 

63 

N 
s 

2800 

10000 

2600 

250 

31 

P  (kips) 

5.5 

15.0 

32.0 

39.8 

45.0 

Lorraine 

Nb 

5200 

18000 

4700 

460 

57 

N 
s 

2600 

9000 

2300 

230 

28 

P  (kips) 

8.2 

18.7 

34.2 

46.2 

75.0 

Aberdeen 

Nb 

4900 

17000 

4400 

430 

54 

N 
s 

2400 

8500 

2200 

220 

27 

P  (kips) 

8.9 

32.8 

38.3 

40.5 

48.0 

FAST 

Nb 

3000 

11000 

2700 

270 

33 

N 
s 

1500 

5300 

1400 

130 

17 

87 


The  equivalent  track  loadings  for  track  sections  having  wood  ties  are 
shown  in  Fig.  29a.   The  mean  dynamic  wheel  loads  were  used  on  three  of  the 
axles.   The  fourth  axle  load  was  based  upon  the  wheel  loads  determined  for 
the  50,  10,  1  and  0.1  percent  exceeded  levels.   The  variable  wheel  loadings 
were  placed  second  in  the  sequence.   The  load  positions  shown  in  Fig.  29a  were 
used  for  the  Leeds  wood  and  FAST  22  wood  tie  sections,  but  the  magnitudes  of 
the  mean  and  variable  dynamic  loads  differed  as  indicated  in  Table  18. 

The  equivalent  track  loadings  for  the  sites  having  concrete  ties  are 
shown  in  Fig.  29b.   Again,  the  load  sequence  used  was  the  same  for  all  con- 
crete tie  sections,  although  the  magnitudes  of  the  applied  dynamic  loads  used 
did  vary  among  the  sites,  according  to  Table  18. 
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CHAPTER  6.   METHODOLOGY  FOR  DEFORMATION  PREDICTIONS 


6.1   Stress  Path  Determination 

The  approach  to  the  prediction  of  the  permanent  deformation  that  accumu- 
lates in  track  due  to  repeated  traffic  loadings  is  based  upon  the  conceptual 
ideas  developed  primarily  in  the  pavement  area.   A  logical  extension  of  the 
pavement  approach  has  been  investigated  for  track  [22,42]   and  will  serve  as 
the  basis  for  the  methodology  presented  in  this  volume. 

The  methodology  for  the  prediction  of  track  permanent  deformations  starts 
by  determination  of  the  stress  states  at  various  locations  within  the  track 
foundation.   These  stress  states  derive  from  two  sources:   the  initial  stres- 
ses due  to  the  unloaded  track  condition,  and  the  incremental  stresses  result- 
ing from  the  imposed  train  loadings.   The  initial  vertical  geostatic  stresses 
are  determined  by  the  weights  of  the  track  superstructure  and  the  soil.   The 
initial  horizontal  geostatic  stresses  are  related  to  the  initial  vertical 
stresses  through  the  K  parameter. 

The  incremental    stresses  are  determined  using  the  GEOTRACK  program. 
The  GEOTRACK  program  computes  the  complete  incremental  stress  tensor  for  var- 
ious locations  in  the  track  formation.   These  incremental  stresses  are  added 
to  the  initial  geostatic  stresses  to  determine  the  final  three-dimensional 
stress  states  at  various  locations.   These  final  three-dimensional  stresses 
are  converted  into  equivalent  triaxial  stress  states  using  the  concept  of 
stress  invariants  [22,47].   The  resulting  triaxial  stress  states  are  then  used 
to  determine  the  representative  stress  paths  for  points  within  the  roadbed. 

A  typical  stress  path  for  a  point  in  the  ballast  directly  below  a  moving 
load,  as  calculated  using  the  invariant  formulation,  is  shown  in  Fig.  30. 
Point  A  represents  the  initial  stress  state.   Since  the  K  value  for  the  bal- 
last layers  is  greater  than  one,  the  initial  shear  stress,  q  .  ,  is  negative. 
Point  B  represents  the  equivalent  stresses  at  a  point  directly  below 

a  wheel  load.   The  cyclic  stress  path  for  the  point  in  the  roadbed  due  to  the 
moving  load  thus  is  given  by  the  line  AB.   Two  alternatives  were  considered 
for  representing  path  AB  as  an  equivalent  constant  confining  pressure  triaxial 
stress  path.   First  is  stress  path  CB,  which  matches  only  the  maximum  shear 

stress,  q    ,  and  the  minimum  or  initial  shear  stress,  q_.  .   Second  is  stress 

,  ~„  max  ,  .  ,      ,  —    •    , ,  win     .      .       , 

path  DE,     which  matches  the  mean  stress,  p  ,  as  well     as  the  maximum  and 

minimum  shear  stresses.   Since  the  actual      traffic  consists  of  a  mix  of 
wheel  loads,  these  stress  paths,  as  shown  in  Fig.  30,  will  be  required  for  a 
variety  of  loading  conditions. 

The  depths  within  the  track  foundation  at  which  the  representative  stress 
paths  should  be  calculated  are  not  uniquely  defined.   The  approach  taken  in 
the  pavement  area  [48]  has  been  to  calculate  the  triaxial  stress  states  at  the 
center  of  the  various  layer  subdivisions  to  determine  the  deformation  contri- 
butions from  each  of  these  layers.   This  may  be  appropriate  when  predictions 
are  to  be  made  for  newly  built  facilities,  but  when  dealing  with  systems  such 
as  maintained  or  rebuilt  track,  the  contributions  of  the  undisturbed  lower 
layers  are  often  negligible  and  can  therefore  be  ignored.   This  will  be  ex- 
plained more  fully  in  a  subsequent  section. 
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Figure  30.  Equivalent  Stress  Paths  for  Roadbed  Points 
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6.2  Strain  Determination 

First  Cycle  Strains 

Once  the  stress  paths  have  been  established  for  the  various  combinations 
of  dynamic  wheel  loadings,  the  next  operation  is  to  determine  the  strains 
that  would  be  expected  to  develop  due  to  repeated  applications  of  the  prescri- 
bed stresses.   This  is  done  using  the  laboratory  results  from  the  various  test 
series  described  in  Chapter  3. 

The  deformations  that  accumulate  in  ballast  due  to  repeated  application 
of  vertical  stresses  were  shown  to  be  related,  in  general,  to  the  permanent 
strain  or  deformation  that  occurred  in  the  first  cycle  of  load.   Thus,  the 
first  step  in  predicting  the  track  deformations  is  to  determine  the  strains 
that  are  expected  to  develop  from  the  first  load  cycle.   Since  the  stress 
paths  for  the  ballast  layers  involve  shear  stress  reversals,  the  accurate  de- 
termination of  these  first  cycle  strains  is  difficult. 

The  first  cycle  strain  for  one-way  repeated  loading  has  been  shown  to  be 
related  to  the  ratio  of  maximum  applied  shear  stress,  q    ,  to  failure  shear 
stress,  qf.   For  tests  involving  shear  stress  reversal,    the  stress  ratio 
used  is  Aq/qf,  where  Aq  is  the  difference  between  the  maximum  and  minimum 
shear  stresses. 

The  relationships  between  first  cycle  strains  and  stress  ratios  for  the 
ballast  materials  are  shown  in  Fig.  31,  both  for  tests  involving  one-way  load- 
ing and  shear  stress  reversals.   The  shear  stress  reversal  data  were  obtained 
from  the  box  tests  and  triaxial  tests.   All  of  the  data  with  a  cyclic  stress 
ratio  greater  than  1.0  are  from  the  box  tests.   The  results  for  all  of  the 
field  site  ballasts  tested,  as  well  as  the  FAST  Wyoming  granite,  were  used  to 
define  the  trends  shown  in  the  figure.   In  each  case,  the  permanent  strain 
after  the  first  cycle,  z    ,  has  been  divided  by  the  effective  confining  pres- 
sure, a_,  in  an  attempt   to  normalize  the  existing  data  so  they  could  be  ex- 
trapolated to  pressures  other  than  those  used  for  the  experiments. 
The  effective  confining  pressures  and  stress  ratios  for  the  box  tests  were  cal- 
culated based  on  the  measured  residual  horizontal  stresses  in  the  box  tests. 

For  each  of  the  field  site  predictions,  the  first  cycle  strains  will  be 
estimated  from  the  curves  shown  in  Fig.  31  for  four  stress  ratios.   These  are: 

Ratio  Stress  Path  Strain  Curve 


qnax/q 


CB  One-Way  Uncompacted 


q    /q^  DE  One-Way  Uncompacted 

max  f 

Aq/q  CB  Uncompacted  Stress  Reversal 

Aq/qf  DE  Uncompacted  Stress  Reversal 

Strain  Accumulation  Rates 

The  rate  of  strain  accumulation  with  repeated  load  cycles  for  the  ballast 
materials  was  previously  discussed  in  Chapter  3.   The  ballast  constants,  C, 
used  for  the  predictions  will  be  both  the  individual  results  given  in  Table  15 
for  one-way  loading  and  the  uncompacted  constant,  C,  ,  derived  from  the  box 
test  results.   These  soil  constants,  together  with    the  permanent  strain  af- 
ter the  first  cycle,  will  be  used  in  predictive  Equation  16  to  determine  the 
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Figure  31.   Normalized  First  Cycle  Strains  as  a  Function  of  Stress  Ratios 
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strain  after  additional  cycles  for  each  of  the  field  site  stress  paths. 
Superposition  of  Strains 

The  cumulative  strains  due  to  the  mix  of  wheel  loads  will  be  handled  on 
the  basis  of  a  cumulative  damage  relationship  similar  to  Miner's  rule.  This 
type  of  cumulative  damage  rule  has  been  frequently  used  in  fatigue  studies 
for  asphalt  pavement  and  pavement  deformation  methodologies  by  a  number  of 
researchers  [49,50,51].   The  validity  of  this  approach  was  established  for 
the  ballast  on  the  basis  of  experimental  results  from  the  staged  testing  pro- 
gram. 

The  basic  procedure  for  application  of  Miner's  rule  starts  with  the  ob- 
servation that  after  N.  cycles  at  stress  level  i,  the  amount  of  permanent 
strain,  e  (i) ,  will     develop.   An  equivalent  number  of  cycles,  N.*,  can  be 

calculated  such  that  the  permanent  strain  after  N.*     J  cycles 
is  equal  to  the  strain  already  developed  after  N.  cycles  at  *•  stress  level  i. 
For  the  ballast  strains,  this  can  be  expressed    as  follows: 

1.  Let  £,  .  and  £..  .  be  the  strain  developed  at  N.  =1  and  N.  =1,  res- 
pectively.  J 

2.  After  N.  cycles  at  stress  level  i,  the  permanent  strain  will  be 
given  by 

eN(i)  =  e1±(l  +  C  log  N±)   ,  (32) 

and  after  N.  cycles  at  stress  level  j,  the  permanent  strain  will  be  given  by 

eN(j)  =  e  (1  +  C  log  N  )   .  (33) 

3.  The  equivalent  number  of  cycles  at  stress  level  j  that  would  have 

been  necessary  to  cause  strains  equal  to  eN(i)  is  found  by  equating  e«(i)  to 

£„,(j)>  in  which  case  N.  becomes  N.*.   The   result  is 
N  J         J 


N.*  =  10 
J 


re*T(i) 


elj 


-  1 


(34) 


4.   The  total  cumulative  strain  after  N.  cycles  of  stress  level  i  plus 
an  additional  N.  applications  of  stress       level  j  is  then  given  by 

e  (i+j)  =  E   [1  +  C  log  (N.*  +  N.)]    .  (35) 

The  method  is  illustrated  in  Fig.  32  for  a  case  with  two  load  levels,  the 
second  of  which  is  greater  than  the  first.   The  procedure  may  be  repeated  for 
all  other  successive  stress  levels,  and  can  even  be  used  for  returning  to 
stress  level  i. 

For  the  field  site  predictions,  the  first  cycle  strains  for  the  50,  10, 
1  and  0.1  percent  exceeded  load  levels  will  be  used,  along  with  the  prescribed 
number  of  cycles  at  each  load  level,  as  given  in  Chapter  5.   However,  for  com- 
parison with  the  field  measurements,  the  calculated  strains  are  plotted  as  a 
function  of  the  actual  number  of  cycles  (or  MGT) ,  not  the  equivalent  number  of 
cycles  used  for  the  superpositions. 
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6.3  Deformation  Determination 

The  calculated  strains  in  each  layer  must  be  multiplied  by  the  layer 
thickness  to  get  the  layer  compression.   All  of  the  layer  compressions  are 
then  summed  to  determine  the  track  vertical  deformation.   However,  the  track 
maintenance  operation  disturbs  only  a  shallow  depth  of  ballast,  at  most  6  to 
8  inches  below  the  base  of  the  ties.   The  layers  and  materials  below  this 
depth  which  have  been  subjected  to  years  of  traffic  and  hundreds  of  MGT 
should  contribute  little  to  the  total  track  settlement,  unless  pumping  or 
embankment  instability  exists. 

The  stresses  used  for  all  of  the  revenue  field  site  predictions  there- 
fore will  be  those  estimated  at  depths  of  3  to  4  inches  below  the  tie,  rep- 
resenting a  layer  thickness  of  6  to  8  inches.   The  contribution  to  settlement 
from  lower  layers  will  be  ignored. 

The  FAST  track  was  rebuilt  and  the  ballast  layer  was  completely  replaced 
in  1979.   Then  a  major  maintenance  operation  was  done  after  about  93  MGT 
traffic.   Thus,  two  sets  of  predictions  are  needed  for  FAST.   The  first  will 
be  for  the  first  93  MGT  traffic  after  rebuild,  considering  deformations  that 
develop  in  the  entire  15-in.  ballast  layer.   The  second  set  of  FAST  predic- 
tions will  be  for  the  traffic  after  the  maintenance  at  93  MGT.   Since  the 
maintenance  affected  at  most  the  top  half  of  the  total  ballast  layer,  this 
prediction  will  consider  only  the  7.5-in.  depth  of  ballast  assumed  to  be 
affected  by  the  disturbance.   The  subgrade  materials  at  FAST  had  been  subjec- 
ted to  over  400  MGT  traffic  prior  to  the  rebuild.   Thus  the  subgrade  contri- 
bution to  settlement  after  the  rebuild  will  be  neglected.   The  errors  result- 
ing from  neglecting  further  subgrade  deformations  should  be  small,  and  cer- 
tainly less  than  the  uncertainty  in  the  method  used  for  the  predictions. 
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CHAPTER  7.   PERMANENT  DEFORMATION  EVALUATION 

7.1  Introduction 

The  track  settlements  at  the  revenue  sites  were  measured  by  BCL  using 
optical  surveys.   For  these  vertical  survey-to-benchmark,  readings,  measurements 
were  taken  at  the  top  of  the  railhead,  at  10-ft  longitudinal  intervals,  at  ad- 
jacent locations  on  each  rail.   At  the  Aberdeen  site,  inductance  strain  coils 
were  also  installed  in  the  ballast  to  determine  the  contribution  of  the  bal- 
last layer  to  the  total  vertical  settlement.   These  coils  were  placed  in  pairs 
beneath  the  rail  seats  under  three  ties.   One  coil  in  each  pair  was  placed  at 
a  depth  of  approximately  13  in.  below  the  bottom  of  the  tie,  in  the  ballast 
layer.   The  other  coil  was  at  the  bottom  of  the  tie.   This  configuration  meas- 
ured the  change  in  upper  ballast  thickness  relative  to  the  initial  thickness 
(about  13  in.)  prior  to  maintenance  when  the  coils  were  installed.   The  princi- 
ple of  operation  for  this  instrumentation  can  be  found  in  Ref.  [52],  along  with 
the  calibration  procedures.   Measurements  of  the  deformations  in  the  track  at 
FAST  also  have  been  recorded  using  vertical  survey  and  foundation  instrumenta- 
tion.  Only  the  measurements  since  the  1979  rebuild  will  be  considered  for  cor- 
relation with  the  revenue  service  track  results.   Measurements  from  the  earlier 
FAST  experiments  prior  to  the  rebuild  have  been  previously  reported  [22,46,53]. 

This  chapter  will  present  the  track  settlement  results  obtained  from  the 
field  sites  and  compare  the  measured  responses  with  the  predictions  made  using 
the  methodology  outlined  in  Chapter  6. 

7.2  Field  Measurements 

Revenue  Sites 

The  correlation  plan  called  for  survey  measurements  to  be  taken  immediate- 
ly after  the  track  surfacing,  prior  to  any  traffic,  followed  by  readings  taken 
at  about  6-month  intervals.   However,  in  some  cases  as  much  as  7  to  10  days 
elapsed,  representing  up  to  0.5  MGT  traffic,  before  the  zero  reference  survey 
was  made.   This  presented  a  significant  uncertainty  in  the  magnitude  of  the 
initial  settlement,  when  relatively  large  displacements  are  likely  to  occur. 

The  actual  railhead  elevations  obtained  by  the  survey  at  each  site  were 
detrended  by  BCL  to  correct  for  the  slight  grades  that  were  present  in  the 
track  vertical  alignments.   Tbese  detrended  profiles  were  referred  to  as  the 
baseline  profiles  for  the  reference  surveys.   For  example,  if  there  was  a  con- 
stant gradient  in  track  elevation  over  the  test  section,  linearly  proportioned 
corrections  to  the  railhead  elevations  were  made  throughout  the  test  section. 
The  initial  corrected  baselines  were  used  as  a  reference  for  the  follow-up 
surveys  to  determine  the  profile  change.   The  mean  detrended  rail  elevations 
from  the  follow-up  surveys  were  subtracted  from  the  mean  initial  rail  eleva- 
tion (baseline)  to  obtain  the  track  settlements. 

The  average  settlement  at  each  of  the  field  sites  as  a  function  of  MGT 
after  maintenance  is  given  in  Table  21.   The  zero  settlement  at  each  site,  as 
previously  explained,  represents  the  average  detrended  reference  elevation  of 
the  test  sections.   The  standard  deviations  of  the  detrended  rail  profiles  for 
the  reference  survey  and  the  follow-up  surveys  are  also  given  in  Table  21. 
This  standard  deviation  can  be  interpreted  as  a  measure  of  track  roughness. 
These  track  roughness  values  did  not  appear  to  be  correlated  to  the  magnitudes 
of  track  settlements. 

97 


4-1 


u    O    4J 


CU 


3 
W 

•J 
CJ 
PQ 

I 


s 

CU 

J-l 

3 
CO 

co 
<u 

S 


d 
cu 

J-i 
•U 
CU 

w 

CJ 

co 

H 


CO 
O 

•H 

4-> 

u 

CU 

> 

TJ 

CU 
U 

o 

CU 

u 

u 
o 
o 

d 


(N 

CU 
tH 

CO 

H 


CO 

d 

ni 

a 

(1 

■H 

43 
Ml 

-i 

■H 

o 

> 

P4 

a) 

O 

ri 

■H 

T3 

rl 

IU 

0 

hfl 

T3 

n 

a 

01 

cfl 

x: 

4J 

u 

on 

o 

^v 

u 

• 

d 

(U 

•H 

> 

vy 

■H 

4-1 

CO 

CO 

CU 

rH 

3 

CO 

iH 

pS 

CO 

> 

en 

4J 

CU 

0 

CJ 

<D 

d 

0 

cu 

OJ 

u 

rH 

CU 

U 

4-1 

■u 

CU 

CU 

pei 

oo 

d 

o 

rH 

•H 

•H 

/— s 

U 

4-1 

CO 

• 

n 

CO 

Pd 

d 

•H 

•H 

ca 

> 

T3 

v^ 

[fl 

CU 

CU 

a) 

Q 

T3 

CO 

C 

d 

CU 

,C 

13 

CU 

rH 

00 

l-i 

n 

•H 

3 

cO 

4-1 

4-4 

O 

•a 

CU 

o 

Pi 

d 

Q 

n 

CO 

CU 

4-1 

4-1 

C/3 

o 

<u 

ad 

a 

C 

c 

•H 

■H 

CJ 

en 

cfl 

4-1 

H 

u 

CJ 

3 

S 

C/J 

cu 
•u 


C/J 


CN 

-a- 

CO 

CO 

CO 

vO 

co 

r* 

CN 

CM 

1 

o 

o 

I 

o 

O 

o 

1 

o 

o 

o 

1 

O 

O 

1 

• 

• 

1 

• 

• 

• 

1 

• 

• 

• 

1 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CT> 

00 

CO 

P>. 

o 

o 

o 

a\ 

CO 

r-» 

vO 

m 

r*-. 

r*. 

CM 

CM 

en 

CN 

CO 

CO 

CO 

iH 

CN 

CM 

CN 

CM 

CM 

CM 

O 

o 

o 

O 

o 

o 

o 

O 

o 

o 

O 

O 

O 

o 

o 

-3- 

v£> 

o 

CX\ 

>* 

r^. 

O 

o 

00 

C^ 

O 

r» 

O 

o 

«tf 

sr 

o 

o 

<fr 

«* 

o 

CO 

iH 

CM 

o 

rH 

m 

o 

O 

o 

o 

o 

o 

o 

o 

o 

r-f 

rH 

o 

O 

© 

■K 

« 

He 

■K 

m 

m 

m 

VO 

iH 

o 

vO 

O 

m 

m 

O 

sr 

o> 

sf 

o 

CN 

CN- 

O 

iH 

<r 

vO 

O 

sr 

vO 

sr 

o 

o 

>3- 

iH 

CN 

CU 

•u 

CU 

u 

iH 

CM 

CO 

CO 

m 

vO 

tH 

CM 

T3 

CJ 

O 

d 

O 

o 

* 

CJ 

CU 

d 

d 

CU 

M 

n 

•H 

CU 

CO 

co 

cO 

T3 

T3 

•o 

u 

U 

CU 

CU 

u 

CU 

CU 

CU 

o 

jq 

J 

-J 

i-3 

< 

CU 

3 

rH 
CO 
> 

CU 
CJ 

d 

CU 

u 

CU 
4-1 
CU 

1-1 

u 
o 


c 

M 


98  - 


The  standard  deviation  of  the  track  settlement  as  referenced  to  the  ini- 
tial baseline  profiles  is  also  given  in  Table  21.   The  standard  deviation  val- 
ues were  computed  from  the  variations  between  the  initial  baseline  and  the  sub- 
sequent profiles,  as  supplied  by  BCL.   These  values  were  small,  indicating 
little  change  in  track  roughness  with  settlement. 

Since  the  zero  reference  surveys  at  the  Leeds  sites  and  at  Aberdeen  were 
taken  after  some  amount  of  traffic  following  maintenance,  a  correction  had  to 
be  estimated  to  account  for  the  missing  initial  settlements  for  these  sites. 
Fortunately,  the  strain  coils  installed  in  the  ballast  at  the  Aberdeen  site 
were  monitored  frequently  during  the  first  field  site  visit.   These  strain 
coils  allowed  a  direct  measurement  of  the  amount  of  raise  at  the  site  and  sub- 
sequent deformations.   However,  the  strain  coil  measurements  would  probably 
give  an  overestimate  of  the  ballast  deformations  throughout  the  entire  track 
section.   The  reason  is  that  a  13-in.  depth  of  ballast  was  excavated  and  re- 
placed when  the  coils  were  installed.   Although  the  ballast  in  the  areas  under 
the  ties  and  in  the  cribs  was  compacted  with  a  portable  mechanical  tamper,  the 
replaced  material  may  not  have  been  as  compact  and  stiff  as  the  undisturbed 
ballast  at  locations  without  strain  coils.   The  effects  of  this  ballast  distur- 
bance at  the  strain  coil  locations  could  have  resulted  in  greater  track  settle- 
ment at  these  locations  than  at  the  remaining  tie  locations. 

The  initial  deformations  that  developed  in  the  upper  13  in.  of  ballast  at 
the  Aberdeen  site  are  shown  in  Fig.  33.   The  coils  indicated  a  surface  raise 
of  about  0.1  in.,  followed  by  a  rapid  initial  settlement.   The  reference  zero 
survey  was  made  10  days  after  the  surfacing.   The  ballast  strain  coils  indica- 
ted that  there  had  been  about  0.20  in.  of  vertical  deformation  in  the  ballast 
since  the  surfacing  due  to  the  first  0.4  MGT  traffic.   Therefore,  the  BCL  sur- 
veys had  to  be  adjusted  by  this  amount  to  correct  for  the  missed  initial  defor- 
mations. 

The  similarity  between  the  trends  of  the  corrected  survey  measurements  and 
the  bal^ajt  strain  measurements  supported  the  conclusion  that  the  subgrade  and 
deeper  Laxlast  settlements  could  be  neglected. 

The  two  Leeds  sites  did  not  have  strain  coils  embedded  in  the  ballast,  so 
the  corrections  to  the  survey  measurements  at  these  sites  were  uncertain.   The 
amount  of  initial  settlement  at  these  sites,  where  the  raise  was  about  1-%   in., 
could  only  be  estimated.   Data  from  FAST  indicated  that  about  .25  in.  of  set- 
tlement occurred  after  0.3  MGT  of  traffic  had  passed  over  section  22  following 
a  comparable  raise.   Based  on  this  information,  the  correction  to  the  Leeds 
survey  measurements  was  estimated  at  0.2  in.   The  actual  survey  data  will  be 
used  for  the  Leeds  sites,  and  the  estimated  corrections  shown  as  an  increase. 
Thus  the  survey  measurements  represent  a  lower  bound  to  the  actual  track  set- 
tlements.  For  the  Leeds  sites,  the  probable  corrections  for  initial  settle- 
ments were  on  the  same  order  as  the  subsequent  settlements  over  the  remaining 
18  months. 

The  survey  was  done  directly  after  maintenance  at  the  Lorraine  site.   How- 
ever, an  interval  of  over  30  MGT  elapsed  from  the  zero  reference  point  to  the 
first  follow-up  survey.   Thus,  the  shape  of  the  settlement  curve  could  not  be 
defined,  even  though  the  end  points  were  known  for  this  traffic  interval. 

FAST  Measurements 

Measurements  of  vertical  track  settlement  at  FAST  section  22  were  obtained 
from  both  embedded  foundation  instrumentation  and  track  surveys.   Details  of 
the  FAST  22  instrumentation  used  after  the  1979  rebuild  can  be  found  in  Ref. 
[54]. 
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The  measured  ballast  and  subgrade  deformations  after  the  rebuild,  as  ob- 
tained from  the  instrumentation,  are  shown  in  Fig.  34a  for  the  wood  tie  sec- 
tion.  The  results  from  the  instrumentation  measurements  for  the  concrete  tie 
section  of  FAST  section  22  are  shown  in  Fig.  34b.   The  vertical  lines  on  these 
figures  represent  ±  one  standard  deviation  of  the  measured  results.   Analyses 
of  the  ballast  and  subgrade  data  shown  in  Fig.  34  indicated  that  there  were  no 
statistically  significant  differences  in  the  measured  ballast  deformations  and 
in  the  measured  subgrade  deformations  between  the  wood  and  concrete  sections 
of  FAST  section  22  after  the  1979  rebuild. 

A  similar  conclusion  regarding  the  ballast  deformations  was  found  from 
the  previous  FAST  static  data  from  sections  17,  18,  and  20.   The  average  bal- 
last deformations  measured  at  FAST  sections  17,  18,  and  20  prior  to  the  rebuild, 
along  with  the  average  ballast  deformations  in  section  22  measured  after  the 
1979  rebuild,  are  compared  in  Fig.  35.   Both  sets  of  measurements  show  an  ini- 
tially high  rate  of  deformation  accumulation.   Considering  the  variability 
associated  with  each  set  of  observations,  it  was  concluded  that  the  deforma- 
tions that  developed  in  the  earlier  FAST  experiments  were  similar  to  those  ob- 
served immediately  after  the  1979  rebuild. 

There  is  some  uncertainty  about  the  ballast  strain  data  taken  after  the 
rebuild.   The  track  rebuild  construction  specifications  indicated  that  the 
ballast  was  to  be  placed  in  2-  to  4-in.  lifts,  followed  by  tamping  and  ten 
passes  of  a  locomotive  with  at  least  10  loaded  100-ton  hopper  cars.   This  pro- 
cedure was  to  be  repeated  until  the  final  grade  was  attained.   Subsequent  in- 
formation regarding  the  rebuild  [6]  reported  6-in.  lifts  followed  by  tamping 
and  regulating,  with  no  mention  of  the  track  loadings.   Since  this  incremental 
construction  and  subsequent  loading  resulted  in  some  ballast  compaction,  as 
was  intended  by  the  rebuild  planners,  even  the  FAST  ballast  strain  coil  data 
from  the  0  to  93  MGT  interval  contain  some  degree  of  uncertainty  in  interpre- 
tation.  The  measured  ballast  strains  shown  in  Fig.  34  represent  lower  bounds 
to  th:  expected  strains  developed  in  a  newly  constructed  track. 

Figure  36  shows  the  cumulative  subgrade  deformations  measured  in  wood  and 
concrete  tie  sections  of  FAST,  since  the  beginning  of  the  FAST  track  operation 
and  after  the  1979  rebuild.   The  wood  tie  section  subgrade  deformations  after 
the  rebuild  were  significantly  larger  than  those  measured  in  the  concrete  tie 
section.   Prior  to  the  1979  rebuild,  the  cumulative  subgrade  deformations  in 
concrete  tie  section  17  of  FAST  were  larger  than  in  the  wood  tie  sections  18 
and  20. 

It  is  unclear  why  the  subgrade  deformations  after  the  1979  rebuild  should 
show  the  significant  changes  indicated  in  Fig.  36.  Only  the  top  few  inches  of 
the  subgrade  apparently  were  disturbed  during  the  rebuild  and  the  remainder  had 
been  subjected  to  over  400  MGT  traffic.  Extrapolation  of  the  trends  from  both 
wood  and  concrete  sections  prior  to  the  rebuild  suggests  that  the  subgrade  de- 
formations due  to  an  additional  100  MGT  after  the  rebuild  should  have  been 
small.   Thus  the  rebuild  subgrade  extensometer  data  appear  to  be  misleading. 

There  are  other  factors  that  could  have  resulted  in  the  presumably  high 
measured  subgrade  deformations  after  the  rebuild.   The  actual  construction 
operations  may  have  disturbed  the  subgrade  to  a  depth  greater  than  the  top  few 
inches.   This  could  have  resulted  from  deeper  undercutting  than  planned,  or 
from  the  vibratory  compaction  used  during  the  final  subgrade  preparation.   Sig- 
nificant vibrations  in  the  sandy  subgrade  may  have  relieved  the  previous  stress 
history  effects,  making  the  material  more  susceptible  to  additional  cyclic-in- 
duced densification.   The  effects  of  vibration  on  the  alteration  of  any  pre- 
vious subgrade  structure  would  be  very  difficult  to  quantitatively  interpret. 
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Figure  34.  Ballast  and  Subgrade  Deformations  from  Instrumentation  at 
FAST  Section  22  after  1979  Rebuild 
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Figure  36.   Measured  Subgrade  Deformations  from  FAST  Instrumented  Sect 
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Instrument  reading  errors,  data  reduction  problems,  or  calibration  errors 
are  also  distinct  possibilities.   All  aspects  of  geotechnical  instrumentation 
require  engineering  judgment  in  the  final  interpretation.   Several  years  of 
experience  directly  related  to  the  FAST  instrumentation  program  have  lead  to 
the  conclusion  that  these  rebuild  subgrade  deformation  measurements  should 
be  disregarded,  or  at  a  minimum,  given  little  emphasis.   The  contributions 
of  the  subgrade  deformation  will  be  disregarded  for  the  prediction  comparisons 
for  FAST. 

Settlement  survey  measurements  were  also  made  at  FAST  section  22  after  the 
1979  rebuild.  These  measurements  (to  be  shown  later)  will  be  used  for  the  com- 
parison with  predicted  track  deformations  for  FAST. 

7.3  Predictions  of  Track  Settlement 

General  Considerations 

The  predictions  of  the  vertical  track  settlements  at  the  revenue  field 
sites  and  FAST  section  22  were  based  on  the  methodology  outlined  in  Chapter  6. 
The  depths  influenced  by  the  revenue  site  maintenance  operations  were  from  6 
to  8  inches  below  the  ties.   The  representative  points  for  these  ballast  depths 
were  thus  taken  as  the  midpoint  depths  for  these  disturbed  ballast  zones. 

Since  the  survey  data  from  FAST  indicated  only  slight  differences  between 
the  wood  and  concrete  track  settlements,  only  predictions  for  the  concrete  tie 
section  of  FAST  section  22  were  made.   However,  these  were  subdivided  into  two 
sets.   The  first  set  was  for  the  0  to  93  MGT  period  and  considered  the  entire 
15-in.  ballast  depth.   The  second  was  for  FAST  section  22  after  the  maintenance 
operations  at  93  MGT.   For  this  set,  a  7.5-in.-deep  disturbed  ballast  layer  was 
assumed. 

The  stress  states  were  estimated  using  GEOTRACK.   The  track  structural 
properties  used  for  the  GEOTRACK  analyses  are  given  in  Table  2.   The  layer  pro- 
files and  divisions  and  subgrade  moduli  for  the  revenue  sites  were  the  same  as 
used  for  the  track  modulus  predictions,  and  are  shown  in  Figs.  22  through  25. 
The  layer  divisions  used  for  FAST  consisted  of  two  7.5-in.  ballast  layers  under- 
lain by  the  silty  sand  subgrade  material.   The  subgrade  moduli  expressions  were 
given  in  Chapter  3  for  the  FAST  section.   The  resilient  strain-shear  stress  bal- 
last formulation,  as  described  in  Chapter  3,  was  used  to  characterize  the  bal- 
last at  all  of  the  sites.   The  multiple-axle  loading  configurations  used  are 
shown  in  Fig.  29,  and  the  appropriate  mean  and  dynamic  wheel  loads  are  given  in 
Tables  17  and  18. 

The  calculated  stress  ratios  and  strains  at  the  first  cycle  for  the  field 
sites  are  given  in  Table  22  for  the  predictions  based  on  stress  path  CB  (Fig. 
30) .   The  calculated  stress  ratios  and  the  first  cycle  strains  for  stress  path 
DE  are  given  in  Table  23.   The  first  cycle  strain  is  obtained  from  the  uncom- 
pacted  data  in  Fig.  31. 

The  conversions  from  MGT  to  number  of  cycles  at  each  of  the  load  levels 
were  explained  in  Chapter  5.   The  method  of  superposition  of  strains  due  to 
the  mixed  loading  conditions  was  explained  in  Chapter  6.   The  ballast  con- 
stants, C,  used  for  each  of  the  predictions  are  given  in  Tables  22  and  23. 

Predictions  based  on  paths  CB  and  DE  using  the  results  from  one-way  re- 
peated loading  will  be  referred  to  as  predictions  CB  (1-way)  and  DE  (1-way) , 
respectively.   Those  predictions  based  on  stress  paths  CB  and  DE  that  used  the 
results  from  the  shear  stress  reversal  and  box  tests  will  be  referred  to  as 
predictions  CB  (SSR)  and  DE  (SSR) ,  respectively. 
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Revenue  Site  Predictions 

The  predicted  track  settlements  for  the  Leeds  wood  tie  section  are  shown 
in  Fig.  37.  The  predictions  made  using  stress  path  DE  gave  the  largest  defor- 
mation and  were  close  to  the  uncorrected  survey  measurements.  The  settlement 
rate  appeared  to  level  off  for  the  second  field  measurement  interval,  and  was 
about  parallel  to  the  predictions.  However,  the  predictions  all  underestima- 
ted the  track  vertical  settlement,  particularly  when  any  adjustments  were  made 
to  correct     the  initial  settlement. 

The  measurements  and  predictions  for  the  Leeds  concrete  tie  track  section 
are  shown  in  Fig.  38.   The  field  survey  measurements  showed  a  very  large  change 
in  track  settlement  from  the  10  to  24  MGT  measurements,  which  was  much  larger 
than  the  initial  change  from  the  0  to  10  MGT  readings.   This  occurred  over  the 
winter  months  and  may  have  been  associated  with  frozen  ballast  or  frost  heave. 
The  two  final  survey  measurements  appear  to  indicate  that  the  track  settlement 
rate  was  leveling  off,  as  was  expected  and  predicted,  but  additional  data  are 
required  to  substantiate  this  conclusion.   The  predicted  values  are  all  much 
lower  than  the  measured  values  after  20  MGT  and  show  a  slower  rate  of  increase 
with  MGT. 

Figure  39  shows  the  measured  survey  track  settlements  and  the  predicted 
track  settlements  for  the  Lorraine  site.   The  field  measurements  indicated  a 
continuous  settlement  of  the  track,  with  no  indication  that  the  deformations 
were  leveling  off  with  increased  MGT.   The  predictions  always  showed  a  lower 
settlement  rate  after  about  10  MGT  traffic.   All  predictions  for  the  Lorraine 
site  overestimated  both  the  initial  settlements  and  the  cumulative  settlements. 
The  first  cycle  deformations  for  the  Lorraine  site  predictions  can  be  deter- 
mined from  Tables  22  and  23,  and  are  not  shown  on  the  figure.   The  first  cycle 
deformations  are  much  smaller  than  the  first  data  points  shown,  since  the 
first  points  plotted  in  the  prediction  figure  represent  about  0.1  MGT. 

The  results  from  the  strain  coil  measurements  and  corrected  survey  meas- 
urements, along  with  the  predictions  of  track  settlement  for  the  Aberdeen  site, 
are  shown  in  Fig.  40.   As  previously  indicated,  the  track  settlement  at  Aber- 
deen apparently  was  due  primarily  to  the  upper  ballast  deformations.   Follow- 
ing a  high  initial  deformation  rate  in  the  first  3  MGT,  the  ballast  deforma- 
tions continued  to  increase  at  a  lower  constant  rate,  showing  no  tendency  to 
level  off  with  increased  traffic  accumulation.   Again,  further  survey  data  or 
strain  coil  measurements  would  have  been  useful  to  help  interpret  longer  term 
trends  past  the  30  MGT  level.   Although  the  measured  values  lie  within  the 
range  of  predicted  values,  the  shapes  of  the  predicted  curves  after  about  5  MGT 
were  quite  different  from  the  observed  settlements. 

FAST  Predictions 

Two  sets  of  predictions  were  made  for  FAST  section  22.   Both  sets  were 
done  for  the  concrete  tie  section  track  parameters.   The  first  set  of  predic- 
tions considered  the  0  to  93  MGT  period  directly  after  the  1979  rebuild.   For 
these  predictions,  the  full  15-in.  depth  of  newly  placed  ballast  was  consider- 
ed.  The  second  set  of  predictions  assumed  that  only  the  top  half  of  the  bal- 
last was  disturbed  due  to  the  maintenance  operation  that  took  place  after  93 
MGT  of  traffic. 

The  measured  and  predicted  vertical  track  settlements  at  FAST  are  shown 
in  Fig.  41.   The  0  to  93  MGT  track  settlements  and  ballast  deformations  both 
showed  an  irregular  increase  from  about  the  30  to  50  MGT  interval,  followed  by 
a  reduced  settlement  rate  up  to  the  time  of  the  surfacing  operations.   The 
shapes  of  both  measured  and  predicted  ballast  deformations  were  similar. 
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However,  the  predictions  for  0  to  93  MGT  were  all  greater  than  the  measured 
ballast  deformations  obtained  from  the  strain  coils.   This  could  be  due  to  the 
compaction  sequences  used  during  the  rebuild  operation,  which  would  have  resul- 
ted in  lower  ballast  strain  under  subsequent  traffic  loading. 

The  difference  between  the  0  to  93  MGT  survey  measurements  and  the  ballast 
deformations  from  the  strain  coil  readings  appeared  consistent  with  the  meas- 
ured subgrade  deformations,  even  though  the  subgrade  deformation  data  are  ques- 
tioned.  This  could  have  resulted  from  compensating  errors.   It  is  known  that 
there  were  some  problems  at  FAST  with  the  0  to  93  MGT  portion  of  the  survey 
data,  with  respect  to  reference  elevations  and  data  reduction  procedures. 

Instrumentation  measurements  of  the  ballast  and  subgrade  cumulative  defor- 
mations were  not  available  after  the  93  MGT  maintenance  at  FAST.   However,  the 
survey  measurements  and  predictions  after  the  surfacing  are  shown  in  Fig.  41. 
Both  the  shapes  of  the  initial  and  longer-term  predicted  settlements  were  simi- 
lar to  the  survey  measurements.   The  magnitudes  were  also  in  general  agreement 
for  the  CB  stress  paths. 

The  measured  track  settlements  from  0  to  93  MGT  after  rebuild  for  the 
strain  coils  and  from  the  surveys  after  the  maintenance  showed  a  similar 
amount  of  initial  track  settlement.   However,  these  data,  as  well  as  the  pre- 
dictions, showed  larger  deformations  from  traffic  after  rebuild  than  were  ob- 
served after  subsequent  maintenance.   This  was  primarily  due  to  the  differ- 
ences in  the  initial  ballast  physical  state  in  the  lower  portion  of  the  ballast 
layers,  even  though  there  was  some  initial  compaction  of  the  ballast  during  the 
rebuild.   This  appears  to  confirm  the  previous  assumptions  that  only  the  upper 
6  to  8  in.  of  the  ballast  layer  should  be  considered  as  disturbed  by  mainte- 
nance.  Considering  the  variability  of  the  FAST  data,  the  agreements  between 
measured  and  predicted  track  settlements  were  good. 

7.4  Discussion 

Stress  Path  J-nsiderations 

The  predictions  made  for  the  field  sites  were  based  on  the  invariant 
transformations  of  the  three-dimensional  stress  tensors  into  axisymmetric,  con- 
stant-conf ining-pressure  triaxial  test  conditions.   The  roadbed  stresses  were 
approximated  by  the  assumed  equivalent  stress  paths  labeled  CB  and  DE  in  Fig. 
30,  since  paths  similar  to  the  one  defined  by  the  moving  load,  path  AB,  could 
not  be  reproduced  using  the  present  laboratory  equipment. 

The  predictions  showed  varying  amounts  of  agreement  with  the  measurements 
that  were  not  clearly  attributable  to  any  of  the  particular  stress  paths.   The 
main  difference  between  the  results  obtained  using  the  two  stress  paths,  CB  or 
DE,  was  in  the  first  cycle  strains  associated  with  these  stress  paths.   The 
choice  of  ballast  constant,  C,  from  either  the  one-way  repeated  loading  tests 
or  from  the  box  test  results  was  of  minor  consequence  compared  to  the  value  of 
first  cycle  strain  used  for  the  predictions.   Considering  the  complex  nature 
of  the  problem  and  the  material  behavior  involved,  and  considering  the  large 
field  variability,  the  magnitude  of  the  predicted  settlements  for  the  revenue 
sites  was  in  reasonable  agreement  with  that  of  the  field  measurements.   Cer- 
tainly the  agreement  was  within  the  range  normally  experienced  in  geotechnical 
predictions.   The  predictions  made  for  the  FAST  site  after  the  93  MGT  mainte- 
nance were  in  particularly  good  agreement  with  the  measured  values,  perhaps 
due  to  the  more  controlled  conditions  present  at  FAST,  especially  in  the  ab- 
sence of  wheel  tread  irregularities. 
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It  would  have  been  possible  to  obtain  closer  agreement  between  the  pre- 
dictions and  the  revenue  site  measurements  by  selecting  other  depth  points  as 
representative  for  the  disturbed  ballast  zones.   This  was  not  attempted  for 
several  reasons.   First,  reasonably  good  agreement  was  obtained  between  meas- 
ured and  predicted  deformations  at  FAST  using  the  midpoints  of  the  two  7.5-in. 
ballast  layers  for  the  0  to  93  MGT  predictions,  followed  by  using  only  the  mid- 
point of  the  upper  layer  as  representing  the  disturbed  ballast  following  main- 
tenance.  Secondly,  the  methodology  used  to  predict  track  settlements  needed 
to  be  applied  consistently  to  all  the  sites  in  order  to  assess  the  results. 
The  use  of  different  representative  depths  in  the  disturbed  ballast  layers  for 
each  site  could  not  be  justified  by  available  information. 

Another  possible  reason  why  the  predicted  deformations  disagreed  with  the 
measurements  is  that  the  actual  stress  paths  for  the  representative  points 
were  not  similar  to  either  stress  paths  AB,  CB  or  DE.   The  box  tests  indicated 
that  the  loaded  and  unloaded  states  resulted  in  essentially  the  same  magnitude 
of  lateral  stress.   In  contrast,  the  GEOTRACK  model  predicts  substantial  in- 
cremental compressive  horizontal  stresses  at  depths  of  4  in.  below  the  tie. 

As  an  example  of  the  differences  that  the  incremental  horizontal  stresses 

have  on  the  stress  paths,  consider  the  calculated  stress  paths  at  =  4  in.  using 

the  50th  percentile  load  from  the  FAST  section  22  analyses,  and  a  K  value  of 

6.   The  initial  vertical  geostatic  stress  at  this  depth  wa£  about  1  psi._  The 

incremental  stresses  from  the  GEOTRACK  calculations  were  Ao  =  12  psi,  Aa  = 

—  xv 

6.8  psi,  and  Aa  =  26  psi.   The  x  and  y  coordinate  directions  Jwere 

parallel  to  the  tie  and  rail,  respectively.   The  z  direction  was  vertical. 
The  calculations  for  the  equivalent  axisymmetric  stresses  based  on  these  incre- 
mental stresses  are  given  in  Table  24.   The  stress  path  equivalencies  assuming 
no  change  in  the  horizontal  stresses  but  the  same  increment  of  vertical  stress 
as  the  box  tests  showed  are  also  given  in  Table  24.   The  incremental  shear 
stresses  on  the  vertical  and  horizontal  planes  were  neglected  for  these  exam- 
ple calculations  since,  in  this  case,  they  were  small. 

The  resulting  stress  paths  are  shown  in  Fig.  42.   The  box  stress  path 
shows  a  peak,  point  at  the  static  failure  envelope,  a  condition  obviously  not 
possible.   If  the  box  test  results  were  indicative  of  the  actual  stress  states, 
it  should  be  clear  from  Fig.  42  that  deformation  predictions  would  be  uncer- 
tain and  very  sensitive  to  the  exact  values  of  stress,  since  the  stresses  cycle 
between  two  failure  states.   This  is  indicated,  in  part,  by  the  permanent 
strain  triaxial  results  in  Chapter  3. 

The  actual  stress  paths  for  the  roadbed  elements  may  lie  somewhere  be- 
tween the  box  test  stress  paths  and  those  predicted  with  the  GEOTRACK  model. 
Static  failure  in  the  conventional  sense  did  not  occur  in  the  box  tests,  due 
to  the  confinement  provided  by  the  rigid  side  and  end  boundary  conditions. 
These  boundary  conditions  were  appropriate  for  measuring  the  K  values  under 
zero  lateral  strain,  but  may  not  have  represented  the  ballast   layer  boundary 
conditions  present  in  the  field.   The  box  tests  did,  however,  reproduce  the 
compaction  effects  of  traffic,  as  was  reflected  by  the  changes  in  ballast  den- 
sity and  plate  bearing  resistances. 

The  predictions  were  not  obtained  for  either  of  the  stress  paths  in  Fig. 
42,  but  for  the  constant-conf ining-pressure  "equivalent"  paths,  as  in  Fig.  30. 
Stress  paths  DE  and  CB  can  be  expected  to  give  different  results  from  those  in 
Fig.  42.   Thus,  if  the  stress  states  actually  lie  somewhere  between  the  GEO- 
TRACK and  box  states  in  Fig.  42,  the  predictions  using  DE  or  CB  should  contain 
some  error. 
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Table  24.  Comparison  of  GEOTRACK  and  Box  Test  Stress  Paths 
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Track  Vibrations 

Another  possible  explanation  why  the  predicted  magnitudes  of  track  settle- 
ment were  generally  less  than  the  measured  values  is  that  the  predictions  and 
laboratory  testing  methods  did  not  consider  any  possible  effects  of  track  vi- 
brations.  Field  measurements  made  by  BCL  at  the  revenue  sites  have  indicated 
that  the  concrete  ties  vibrate  considerably  at  three  principal  bending  modes: 
120,  350  and  650  Hertz,  particularly  in  the  presence  of  wheel  impacts.   The 
actual  vibrations  transmitted  to  the  ballast  were  not  measured,  but  the  coup- 
ling between  the  tie  and  the  top  ballast  particles  is  assumed  to  be  poor  at 
these  frequencies.   Field  measurements  of  the  vibration  frequencies  in  rail- 
road track  foundations  [55]  have  shown  transient  values  on  the  order  of  30  to 
40  Hz  in  the  supporting  soil  layers.   The  conclusion  was  drawn  by  Morgan  [56] 
that  traffic-induced  vibrations  in  the  ballast  could  lead  to  more  rapid  accu- 
mulation of  permanent  deformation  than  would  result  from  repeated  loadings 
alone. 

The  magnitudes  of  the  ballast  deformations  that  could  have  been  attributed 
to  track  vibrations  is  difficult  to  quantify.   It  is  likely  that  vibrations 
coupled  with  high  loads  from  wheel  irregularities  would  be  most  damaging.   Be- 
cause the  vibrations  occur  simultaneously  with,  and  as  a  result  of,  the  repeat- 
ed loadings,  separation  of  the  settlements  due  only  to  vibration  from  those 
due  primarily  to  the  mixed  repeated  loadings  would  be  very  complicated. 

The  Aberdeen  site  could  have  been  subjected  to  vibrations  sufficient  to 
cause  the  difference  between  predicted  and  measured  trends  of  settlement  with 
MGT.   Here  there  were  high-speed  passenger  cars  that  were  shown  to  cause  the 
high  dynamic  loads  as  a  result  of  wheel  tread  irregularities.   Vibrations  trans- 
mitted to  the  ballast  could  easily  cause  deformations  in  excess  of  those  re- 
sulting from  cyclic  loading  alone. 

Although  the  predicted  Lorraine  deformations  were  larger  than  the  meas- 
ured values,  the  shapes  of  the  predicted  and  measured  curves  were  also  differ- 
ent.  The  vibration  effects  at  the  Lorraine  site  were  not  expected  to  be  as  se- 
vere as  Aberdeen.   Without  actual  field  measurements,  however,  differences  in 
the  effects  of  vibration  at  these  sites  is  uncertain.   The  ballast  at  the  Lor- 
raine site  was  heavily  fouled,  but  still  initially  loose  due  to  the  high  raise 
applied  to  the  track  during  the  surfacing.   Clean  ballast  would  be  more  suscep- 
tible to  vibratory  densif ication,  but  laboratory  experiments  that  were  conduct- 
ed indicated  that  fouled,  compacted  ballast  also  can  be  densified  by  vibration 
frequencies  and  accelerations  on  the  order  of  those  anticipated  in  track. 

The  relatively  good  agreement  between  measured  and  predicted  settlements 
for  FAST  are  probably  a  result  of  the  absence  of  significant  vibrations  at 
FAST.   The  consists,  speeds,  and  wheel  conditions  at  FAST  were  very  well  con- 
trolled, compared  to  the  revenue  sites. 

Other  Factors 

In  addition  to  the  ballast  deformations  that  resulted  from  cyclic  load- 
ings and/or  vibration,  there  are  numerous  factors  that  have  not  been  accounted 
for  in  the  prediction  methodology.   These  include  ballast  recompaction,  bal- 
last degradation,  and  subgrade  effects,  as  well  as  possible  seasonal  effects. 

Ballast  Recompaction.   Laboratory  triaxial  tests  did  not  show  significant 
amounts  of  volume  change  under  repeated  loadings.   Thus,  these  tests  did  not 
reproduce  the  ballast  densif ication  or  compaction  measured  in  the  field.   How- 
ever, the  box  tests  did  reproduce  the  compaction  changes  in  the  ballast. 

The  measured  changes  in  ballast  density  due  to  the  maintenance  operations 
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could  be  interpreted  in  a  form  analogous  to  the  comparable  strains  or  volume- 
tric reductions  resulting  from  one-dimensional  compressions.   It  can  be  shown 
that  the  relative  increase  in  density  of  the  ballast  under  zero  lateral  strain 
conditions  is  given  by 


Tj  ~  Yj-   AY 


b 


d    'do    'd      v 


(36) 


Yj     Yj    1  -  e 

do      do        v 

where  Y-,  is  the  dry  density  at  any  time, 
d 

Y,   is  the  initial  dry  density  after  maintenance,  and 

e   is  the  vertical  strain,  assuming  that  densif ication  occurs  with  no 
horizontal  strain. 
From  pre-  and  post-maintenance  ballast  densities,  particularly  under  the  ties, 
the  expected  track  settlements  due  to  recompaction  of  the  ballast  could  be  es- 
timated using  Eq.  36.  The  rate  and  time  history  of  the  expected  settlements 
could  not  be  determined,  since  the  ballast  densities  at  only  the  two  points 
were  known.   Using  the  results  from  the  pre-  and  post-maintenance  BDT  results 
[1,42],  the  vertical  settlements  from  ballast  recompaction  were  estimated. 
The  Leeds  concrete  site  pre-maintenance  densities  were  not  available,  as  pre- 
viously explained.   The  Aberdeen  under-tie  BDT  showed  negligible  changes  in 
ballast  density  due  to  the  maintenance  operations.   Thus,  the  Leeds  wood  sec- 
tion and  Lorraine  results  for  under-tie  tests  will  be  used  to  conceptually 
illustrate  the  settlements  expected  due  to  recompaction. 

The  measured  density  change,  the  disturbed  zone  thickness  and  the  result- 
ing layer  compression  for  the  two  sites  are  as  follows: 


Site 

AYd 
Ydo 

0.094 
0.115 

e 

V 

h 
(in.) 

6 

8 

e  •  h 

V 

(in.) 

>eeds,  Wood 
Lorraine 

0.086 
0.103 

0.52 
0.83 

The  assumed  influence  depths,  h,  of  the  recompacted  zones  were  the  same  as  the 
disturbed  ballast  depths  used  for  the  strain  predictions. 

As  can  be  seen  from  the  measured  track  settlements  shown  in  Figs.  37  and 
39,  the  above  values  of  settlement  due  to  recompaction  of  the  ballasts  are  of 
the  same  order  of  magnitude  as  the  field  settlement  measurements. 

The  mechanisms  of  ballast  recompaction  are  the  same  as  those  affecting  all 
aspects  of  track  settlement.   All  of  these  mechanisms  occur  simultaneously, 
i.e.,  cyclic  loading,  vibration,  and  possible  horizontal  shearing  at  the  tie- 
ballast  interface.   Individual  contributing  factors  cannot  be  separated  easily, 
then  analyzed, and  finally  recombined  in  an  additive  manner. 

Ballast  Degradation.   Breakdown  of  the  ballast  particles  will  cause  settle- 
ment because  of  a  decrease  in  volume  of  the  assembly  of  particles,  similar  to 
the  effect  of  compaction.   The  causes  of  long-term  ballast  degradation  can  de- 
rive from  several  sources.   Particle  breakage  from  the  high  intergranular  con- 
tact stresses  is  one  source.   The  susceptibility  of  different  ballasts  to  such 
breakdown  is  currently  not  well  defined,  nor  are  the  effects  that  such  degra- 
dation would  have  on  track  performance.   Research  is  currently  underway  at 
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UMass  to  investigate  the  breakdown  potential  for  various  ballasts,  including 
some  of  the  materials  used  at  both  the  revenue  sites  and  FAST,  so  that  better 
performance  evaluations  of  ballast  can  be  made. 

Ballast  degradation  also  can  be  affected  by  factors  not  associated  with 
trattic  loadings,  such  as  chemical  attack  and  weathering.  These  factors  af- 
fect different  ballasts  in  different  ways. 

Freeze- thaw  cycles,  in  areas  prone  to  such  seasonal  variations,  could 
also  contribute  to  particle  breakdown.   Ice  intrusion  into  the  small  fissures 
in  individual  ballast  particles  could  lead  to  increased  cracking  and  reduced 
particle  strength.   A  completely  frozen  ballast  condition,  where  the  indivi- 
dual particles  were  "cemented,"  could  lead  to  higher  stresses  in  the  formation 
and  hence  increased  particle  breakage. 

Subgrade  Effects.   The  contribution  of  subgrade  to  the  total  track  settle- 
ment was  neglected  in  the  predictions.   The  relationship  between  the  strain 
coil  measurements  at  Aberdeen  and  the  agreement  between  the  FAST  track  survev 
and  predictions  partially  justified  this  simplification. 

Another  reason  why  the  subgrade  contributions  were  not  included  was  that 
the  measured  permanent  strains  from  the  laboratory  tests  done  on  the  revenue 
site  subgrade  samples  were  inconsistent  with  either  the  measured  or  expected 
settlements.   In  some  cases,  the  strains  developed  in  subgrade  samples  under 
assumed  representative  triaxial  stress  conditions  became  very  large  after  a 
relatively  small  number  of  load  cycles.   This  may  have  been  the  result  of  in- 
correct representative  stress  paths,  sample  disturbance,  the  fact  that  the 
cyclic  loading  conditions  were  undrained,  or  a  combination  of  all  these 
ractors. 

The  undrained  loading  condition  is  reasonable  for  representing  the  train 
loading  effects,  but  continued  cycling  without  allowing  the  dissipation  of  any 
excess  pore  pressure  that  might  develop  is  not  representative  of  the  field  con- 
ditions  The  trains  do  not  pass  continuously  over  the  sites,  but  are  spaced  at 
relatively  long  intervals.   This  allows  pore  pressure  dissipation  to  occur  be- 
tween load- ngs,  and  hence  decreasesthe  resulting  permanent  deformations.   The 
radial  strain  developments  in  the  triaxial  tests  with  constant  lateral  stress 
would  also  lead  to  larger  axial  deformations  than  anticipated  in  the  field 
where  the  subgrade  soils  are  greatly  confined  by  the  surrounding  ground. 

The  extrapolation  of  the  triaxial  test  permanent  strain  data  would  have 
provided  excessively  large  predicted  values  of  subgrade  deformation.   This  was 
particularly  true  for  the  Aberdeen  site,  where  it  was  found  that  there  were 
negligible  subgrade  contributions  to  the  track  settlement.   The  effects  of  par- 
tial or  intermittent  drainage,  plane  strain  versus  triaxial  testing,  and  stress 
path  effects  on  the  subgrade  deformation  characteristics  are  areas  that  require 
further  study.  H 

Other  subgrade  conditions  can,  however,  have  a  large  influence  on  overall 
track  performance.   Factors  such  as  subgrade  consolidation  under  static  embank- 
ment loads  and  subgrade  movements  from  frost  and  moisture  changes  would  be  im- 
portant in  regions  having  subsurface  conditions  susceptible  to  these  mecha- 
nisms.  Even  after  many  hundreds  of  MGT,  soil  embankments  may  continue  to  set- 
tle by  lateral  creep  induced  by  the  repeated  traffic  loads,  even  though  con- 
fined subgrades  are  stable.   The  Lorraine  site,  in  particular,  may  be  suscep- 
tible to  this  condition. 

For  newly  built  track,  subgrade  compaction  from  the  cyclic  traffic  loads 
can  be  significant.  The  increases  in  penetration  resistance  in  the  FAST  sub- 
grade  resulting  from  the  first  400  MGT  clearly  showed  densif ication 


120 


Subgrade  penetration  into  the  ballast  when  inadequate  soil  filters  or  no 
filters  are  used  to  prevent  particle  migration  is  clearly  another  source  of 
settlement.   However,  the  field  sites  involved  in  this  study  did  not  appear 
to  involve  this  problem. 
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CHAPTER  8.   SUMMARY  AND  CONCLUSIONS 

The  objective  of  this  portion  of  the  concrete  crosstie  correlation 
study  was  the  prediction  of  vertical  track  settlement  at  four  field  sites 
Since  the  track  systems  selected  for  analysis  were  in  different  geographical 
locations  and  subjected  to  different  loading  and  environmental  conditions,  it 
ZH"!?  to  perform  a  complete  ballast  physical  state  assessment  and  sub- 
grade  investigation  of  the  track  foundations  present  at  each  site 

h«l  W%bal-?S\Phy8iCai  ftatC8  W6re  identified  "sing  the  plate  load  test, 

till   of  r T  Y   f '     lateral  tle  PUSh  teSt'   A  ^i-tenance  cycle,  consis- 
ting of  raising  and  tamping  the  ballast  layer,  was  performed  at  each  of  the 

ermd  the  resulting  changes  in  ballast  conditions  were  measured.   One  of 
the  primary  effects  of  the  track  maintenance  was  a  loosening  and  disturbance 

IL   ^J,rTTt\C°mVTed   ^  Stable  ballaSt  layer'   The  ballast  structure 
was  disturbed  both  in  the  cribs  and  under  the  ties.   The  amount  of  disturbance 

Zcrlft  ^   tb1.am°Unt  °f  the  raise«   ^  was  found  that  for  trackbeds  sub- 

jected to  a  significant  raise,  say  1.5  to  2  in.,  the  physical  states  of  the 
ballast  immediately  after  the  raise  were  approximately  independent  of  ballast 
type  and  initial  condition.  " 

_    Samples  of  the  ballast  materials  were  recovered  from  the  sites  and  sub- 

from6FA^  rePeate"  i°ad  triaXial  teStlng-   Additionally,  a  granite  ballast 
from  FAST  was  tested  in  triaxial  repeated  load  conditions  with  and  without 
shear  stress  reversals.   Partial  unloading  and  sequential  application  of  re- 
peated loads  were  also  considered.   It  was  found  that  the  resilient  moduli 
values  for  full  unloading  of  compressive  vertical  stresses  could  be  correlated 

for^rXimUm  bU  k  Stref  °n  the  Sample'   Elther  the  lo8-lo8  or  an  arithmetic 
formulation  could  be  used  to  relate  the  compressive  resilient  Young's  modulus 
to  the  maximum  bulk  stress.   For  partial  unloadings,  however,  the  resilient 
modulus  was  ^igher  than  that  measured- under  full  unloading,  and  was  not  unique- 
ly related  to  the  maximum  bulk  stress  on  the  sample. 

-»,  -  luB   reSU^S  from  the  shear  ^ress  reversal  repeated  load  testing  showed 
that  the  resilient  moduli  from  extension  and  compression  stress  states  were 
quite  different.   The  extensional  resilient  moduli  values  were  smaller  than 
the  congressional  moduli  values  under  comparable  stress  conditions.   When  the 
vertical  stresses  were  cycled  between  compression  and  extension  states,  the 
total  resilient  moduli  values  fell  between  the  extensional  and  compressive 
values.   This  is  believed  to  more  closely  represent  actual  track  conditions 
than  prior  techniques,  which  did  not  include  shear  stress  reversals 

A  comparison  of  the  resilient  moduli  of  all  the  field  site  ballasts  also 
leads  to  the  conclusion  that  the  moduli  values  were  not  significantly  depen- 
dent upon  ballast  type. 

The  inelastic  behavior  of  the  ballast  was  investigated  using  the  repeated 
load  triaxial  tests  and  a  box-type  testing  device.   The  permanent  strain  beha- 
vior of  single-staged,  constant-amplitude  tests  without  shear  stress  reversal 
could  be  defined  by  the  first  cycle  strain,  e   and  the  material  constant,  C. 
This  constant  was  independent  of  stress        conditions  for  a  particular  bal- 
last, but  was  not  the  same  for  all  ballasts. 

Staged  tests  on  the  ballast  were  done  involving  sequential  application  of 
variable  magnitude  repeated  loads.   For  these  tests,  the  total  permanent  strain 
developed  in  a  ballast  sample  was  controlled  by  the  maximum  past  load  and  was 
independent  of  the  loading  sequence,  so  long  as  the  number  of  load  repetitions 
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at  each  stress  level  was  the  same  for  each  of  the  tests.   This  means  that 
the  heaviest  wheel  loads,  generally  produced  by  impacts  from  wheel  tread 
irregularities,  cause  most  of  the  settlement.   Subsequent  applications  of  re- 
duced load  as  well  as  partial  unloadings  resulted  in  negligible  increases  in 
permanent  strain. 

When  the  repeated  loads  in  the  triaxial  tests  were  cycled  between  exten- 
sion and  compression,  the  permanent  strains  developed  at  rates  significantly 
greater  than  the  rates  for  one-way  loading.   The  increased  permanent  strain 
development  was  dependent  upon  both  the  compressive  failure  ratio  and  the  ra- 
tio of  cyclic  shear  stress  to  the  maximum  shear  stress  resulting  from  compres- 
sion loading. 

The  permanent  deformation  accumulation  in  the  ballast  box  tests,  which 
were  intended  to  simulate  field  behavior,  was  a  logarithmic  function  of  the 
number  of  load  cycles,  the  compaction  state,  and  the  first  cycle  deforma- 
tions.  The  trends  in  inelastic  behavior  were  similar  to  those  in  the  triax- 
ial tests.   The  relationship  between  permanent  deformation  and  the  log  of 
the  number  of  cycles  was  similar  to  the  trend  of  permanent  deformation  ob- 
served at  FAST. 

The  box  tests  also  showed  that  relatively  large  horizontal  residual 
stresses  can  build  up  in  the  ballast  as  a  result  of  cyclic  loadings.   These 
residual  stresses  can  develop  not  only  under  the  tie,  but  also  in  the  crib 
areas.   This  may  be  a  mechanism  in  the  phenomenon  of  tie  migration  or  skew- 
ing.  The  horizontal  stress  measurements  made  in  the  ballast  box  suggested 
that  the  ballast  layer  in  track  is  subjected  to  extreme  variations  in  cyclic 
stresses.   The  unloaded  stress  states  may  be  close  to  the  limiting  passive 
failure  stress  conditions. 

Measurements  of  vertical  track  modulus  were  made  before  and  after  the 
revenue  site  maintenance  operations.   The  major  factors  contributing  to  the 
magnitude  of  vertical  track  modulus  were  ballast  depth  and  subgrade  stiff- 
ness.  Track  modulus  measurements  made  in  the  track  sections  at  Leeds  having 
wood  txes  spaced  at  19.5  in.  and  concrete  ties  spaced  at  24  in.  were  not  sig- 
nificantly different.   The  field  measurements  did  not  indicate  significant 
changes  in  the  magnitude  of  track  modulus  due  to  maintenance.   Since  the  bal- 
last is  always  stiff  relative  to  the  subgrade  and  since  the  relative  contri- 
bution of  the  ballast  deformation  to  the  overall  elastic  track  deflection  is 
small,  this  result  may  be  expected. 

The  field  measurements  did  show,  however,  that  the  uniformity  of  the  sup- 
port along  the  track  was  not  improved  as  a  result  of  the  surfacings.   Although 
the  property  of  the  ballast  under  individual  ties  was  made  more  uniform  because  of 
the  raising  and  tamping,  as  reflected  by  the  plate  load  tests,  the  uniformity 
of  contact  between  the  ties  and  the  ballast  surface  may  not  be  improved. 

The  GEOTRACK  model,  with  the  stress-dependent  ballast  modulus  formula- 
tion and  constant  subgrade  moduli  determined  for  each  of  the  sites,  was  used 
to  predict  track  deflections  and  track  moduli.   The  6-  to  30-kip  wheel  load 
intervals  were  used  for  the  comparisons  since  the  field  measurements  showed 
some  initial  slack  in  the  tracks.   The  model  predictions  of  vertical  track 
deflection  agreed  well  with,  although  they  tended  to  underestimate,  the  field 
measurements. 

The  original  use  of  track  modulus  was  to  enable  track  designers  and  en- 
gineers to  determine  the  load-carrying  capacity  of  track  in  a  reliable  manner. 
With  structural  components  such  as  concrete  crossties  and  flexible  tie  pads  of 
selectable  stiffness,  the  value  of  the  calculated  track  modulus  may  not  be 
useful  as  a  performance  indicator  if  the  assessment  is  based  only  on  the 
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average  magnitude  of  track  modulus.   The  variability  of  track,  modulus  be- 
tween tie  locations  would  be  a  more  useful  standard  of  track  quality,  since 
the  variability  would  be  a  direct  measure  of  the  uniformity  of  the  track  sup- 
port, which  is  very  important. 

The  measured  dynamic  track  loadings  were  used  to  characterize  the  mixed 
loading  conditions  for  each  of  the  field  sites.   The  dynamic  wheel  load  meas- 
urements and  static  wheel  loads  obtained  from  consist  data  showed  that  the 
mean  static  and  dynamic  loads  were  approximately  equal  for  a  given  field  site. 
The  distribution  of  these  loads  was,  however,  quite  different.   The  dynamic 
loads  below  the  10  percent  level  of  exceedance  were  not  significantly  differ- 
ent from  the  static  loads.   However,  above  this  level,  the  dynamic  loads  were 
much  greater  than  the  static  loads  for  a  given  percent  exceedance  due  to  the 
impacts  resulting  from  wheel  tread  irregularities. 

At  the  0.1  percent  exceedance  levels,  the  dynamic  wheel  loads  at  Leeds, 
Lorraine, and  FAST  were  similar,  averaging  45.5  kips,  with  a  maximum  differ- 
ence of  only  4.5  kips.   However,  the  0.1  percent  exceedance  load  at  Aberdeen 
was  75  kips.   The  very  high  dynamic  loads  measured  at  the  Aberdeen  section 
were  3  to  5  times  greater  than  the  mean  dynamic  load.   These  high  wheel  loads 
are  the  result  of  a  combined  effect  of  wheel  tread  irregularities  and  the 
high  speeds  of  passenger  trains. 

The  loading  environment  at  FAST  was  much  more  uniform  that  at  the  revenue 
sites,  because  of  the  controlled  nature  of  the  FAST  experiments.   While  both  the 
mean  dynamic  and  static  loads  at  FAST  were  on  the  order  of  1.5  to  2  times  the 
average  values  measured  at  the  revenue  sites,  the  low  probability,  large 
loads  were  0.5  to  0.8  times  the  comparable  loads  at  the  revenue  sites. 

A  method  was  developed  to  convert  MGT  into  an  equivalent  number  of  load 
cycles  for  the  ballast  and  subgrade.   The  equivalent  number  of  load  cycles 
per  MGT  for  the  ballast  ranged  from  27,000  to  31,000  at  the  revenue  sites  and 
was  about  17,000  at  FAST.   For  the  subgrade,  the  corresponding  numbers  of  cy- 
cles per  MGT  are  one-half  the  numbers  in  the  ballast. 

Fielc  measurements  of  vertical  track  settlement  were  made  by  BCL  at  each 
of  the  revenue  sites.   These  measurements  were  made  using  optical  survey  tech- 
niques.  The  survey  data  from  the  Aberdeen  site  were  supplemented  by  inductance- 
type  strain  coils  installed  in  the  ballast  layer.   Measurements  of  track  set- 
tlement at  FAST  were  obtained  from  both  optical  surveys  and  embedded  instru- 
mentation.  Estimates  of  initial  settlement  at  the  Leeds  wood  and  concrete  tie 
sections  were  necessary  to  account  for  the  effects  of  about  0.5  MGT  of  traffic 
having  passed  over  the  newly  surfaced  test  sections  prior  to  the  initial  base- 
line surveys.   The  ballast  strain  coils  at  the  Aberdeen  site  allowed  an  accu- 
rate correction  to  be  applied  to  the  baseline  survey  at  that  site. 

Permanent  ballast  strain  data  acquired  using  the  inductance  coils  indica- 
ted that  the  total  track  vertical  deformations  at  Aberdeen  could  be  accounted 
for  by  the  upper  ballast  deformation  alone,  even  though  only  minimal  distur- 
bance of  the  ballast  occurred  there.   This  leads  to  the  conclusion  that  the 
zone  of  ballast  that  was  disturbed  by  the  maintenance  operation  was  the  pri- 
mary source  of  track  settlement,  and  that  the  subgrade  contribution  could  be 
neglected. 

The  settlements  at  the  Lorraine  and  Aberdeen  sites  continued  to  develop 
at  what  appeared  to  be  a  steady  rate  throughout  the  period  of  the  measurements, 
after  the  initial  traffic  following  maintenance.   Although  the  total  vertical 
settlements  at  all  of  the  revenue  sites  continued  to  increase  as  MGT  accumula- 
ted, the  roughness  of  the  track  profiles  did  not  increase. 
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The  ballast  deformations  at  FAST  were  found  to  be  similar  for  both  wood 
and  concrete  tie  sections.   This  was  true  for  the  ballast  measurements  obtain- 
ed during  the  initial  FAST  experiments  and  for  the  ballast  deformations  fol- 
lowing the  1979  rebuild.   The  measurements  of  total  track  settlement  and  bal- 
last layer  deformations  in  FAST  section  22  tended  to  support  the  conclusion 
that  the  subgrade  contributions  to  the  track  settlement  could  be  neglected. 

The  methodology  to  predict  the  inelastic  or  permanent  settlement  of  the 
field  track  sections  used  GEOTRACK  nonlinear  stress  analyses  of  the  field 
site  track  sections  with  stress-dependent  ballast  properties,  representative 
subgrade  moduli,  and  multiple-axle  combinations  of  the  measured  dynamic  wheel 
loads.   The  predicted  stresses  from  the  computer  analyses  for  each  particular 
field  site  and  loading  condition  were  used,  along  with  the  permanent  deforma- 
tion behavior  of  the  ballasts  from  laboratory  tests  to  estimate  the  resulting 
ballast  strains. 

The  mixed  dynamic  field  loading  conditions  were  accounted  for  by  a  formu- 
lation similar  to  Miner's  rule  for  fatigue  analysis.   The  type  of  strain  su- 
perposition was  justified  based  upon  the  results  of  the  staged-test  laboratory 
experiments.   The  settlement  caused  by  the  ballast  was  then  determined  as  the 
product  of  the  ballast  strain  and  the  depth  of  ballast  disturbed  by  mainte- 
nance.  The  disturbed  ballast  depth  was  taken  to  be  only  6  to  8  in.  below  the 
bottom  of  the  tie  for  the  maintained  track  sections.   The  contributions  of 
the  subgrade  to  the  total  track  settlement  were  neglected. 

Although  the  general  magnitudes  of  the  predicted  deformations  agreed 
with  the  measured  values,  the  shapes  of  the  deformation  versus  MGT  predicted 
results  did  not  agree  as  well.   This  was  particularly  true  for  the  revenue 
sites.   Considering  all  of  the  sites,  the  settlement  predictions  made  for 
FAST  section  22  agreed  best  with  the  field  measurements.   This  was  probably 
a  result  of  the  more  controlled  track  loading  and  test  conditions  present  at 
FAST. 

The  method  used  to  determine  the  ballast  strains  incorporated  two  alter- 
native types  of  stress  paths,  one  that  matched  the  peak  shear  stresses  and 
one  that  matched  both  the  peak  shear  stresses  and  mean  stress  of  a  typical 
roaobed  point.   There  were  no  distinct  trends  to  indicate  that  either  of  the 
assumed  stress  paths  based  on  the  GEOTRACK  calculations  led  to  more  consis- 
tent agreement  between  measured  and  predicted  response. 

Uncertainties  in  both  the  first  cycle  strains  and  the  representative 
stress  paths  could  account  for  some  of  the  differences  between  measured  and 
predicted  settlements.   The  results  from  the  box  test  experiments  indicated 
very  small  or  negligible  horizontal  stress  increases  from  load  application, 
even  in  the  top  few  inches  of  the  ballast  layer.   Stress  paths  derived  from 
these  conditions  would  predict  failure  stress  states  in  the  ballast  layer. 
In  contrast,  the  high  K  values  assigned  to  the  ballast  and  the  substantial 
increases  in  horizontal  stresses  in  the  upper  ballast  predicted  by  the  GEO- 
TRACK model  resulted  in  very  different  stress  paths  and  non-failure  condi- 
tions.  Such  uncertainties  in  the  actual  ballast  stress  path  make  accurate 
determination  of  the  first  cycle  strains  difficult.   The  actual  representa- 
tive stress  paths  for  the  track  ballast  layers  may  be  intermediate  to  the 
bounds  created  using  the  box  test  and  GEOTRACK  results. 

Deviations  between  measured  and  predicted  track  settlements  could  have 
resulted  from  other  factors  not  considered  in  the  methodology.   One  of  the 
most  important  factors  not  included  is  thought  to  be  settlement  induced  by  vi- 
brations transmitted  through  the  track  superstructure  to  the  ballast,  as  well 
as  the  extremely  low  probability,  large  loads  encountered  in  revenue  track. 
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Other  factors  that  can  have  varying  degrees  of  influence  on  track  settlement 
and  performance  include  ballast  degradation,  environmental  effects  and  sub- 
grade  responses.   Although  these  last  three  factors  did  not  appear  signifi- 
cant for  the  monitoring  periods  at  the  field  sites  investigated  during  this 
study,  they  could  well  be  significant  at  other  sites. 
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